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ABSTRACT 


That  portion  of  the  line  absorption  eoefficient  for  heated  air  which  depends  only  on  Ihe 
molecular  band  systems ,  i.e. «  not  on  population  numberSi  has  been  calculated  for  each 
line  in  the  frequency  region  of  interest^  with  a  total  of  over  150^000  lines  tabidated  on 
magnetic  tape^ 

With  the  results  for  the  0^  Schumann^ Ruage  system^  a  preliminary  transport  calcu¬ 
lation  has  been  completed  for  a  slab  of  molecidar  os^gen.  The  effect  of  line  widths  on 
transmission  is  illustrated »  and  average  transmissions  are  given  as  a  functimi  of  fre¬ 
quency.  The  present  results  are  compared  with  those  of  a  previous  calculation. 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


( 

1 

1 

I 

I 


II 

I 
1 


2-57-62-2 


CONTENTS 


Sectim  Page 

- EOREWORIX - - - iii- 

ABSTRACT  v 

illustrations  ix 

1  INTROPUCTION  1-1 

2  RADIATION  TRANSPORT  2-1 

2.1  The  Equation  of  Radiation  Transfer  2-1 

2.  2  The  Planck  and  Rosseland  Mean  Absorption  CoefficiintB  2-2 

2. 3  Importance  of  Line  Widths  in  Radiation  Transport  2-5 

3  LINE  ABSORPTION  COEFf  ICIENT  OF  HEATED  AIR  3-1 

3. 1  Expression  for  the  Absorption  Coefficient  of  a  Single  Line  3-1 

3.2  Overlapping  and  Composite  Lines  3-7 

4  THE  LINE  FREQUENCllS  AND  LmE  INTENSITIES  FOR  AIR  4-1 

4. 1  Spectroscopic  Quantum  Numbers  4-1 

4.2  Symmetry  Properties  of  the  Eigenfunctions  4-2 

4. 3  Coupling  of  the  AngiHar  Momenta  4-6 

4.4  Selection  Rules  and  Branches  4-9 

4.5  Intensity  Factors  4-11 

4.6  The  B  -  X  (Schumann>'Runge)  System  of  0„  4-12 

U  ,  &  '  " 

4.7  The  B  -  n  -  A  (First  Positive)  System  of  N2  4-14 

4.8  The  N2  C  ^Hu  -  B  ^Hg  (Second  Positive)  System  4-17 

4. 9  The  B  'Ey  =  X  ^Zg  (First  Negative)  System  of  N2  4-22 

4. 10  The  B  -n  =  X  -H  (Beta)  System  of  NO  4-24 

4. 11  The  A  %  -  X  (Gamma)  System  of  NO  4-28 

5  THE  SACHA  CODES  5=1 

5.1  Magnetic -Tape  Line  Atlases  5-1 

5.2  Average  Transmission  Calculation  5-2 

5.3  Comparison  of  Results  5-3 

5.4  A  Calculation  for  Heated  Air  5-8 

6  REFERENCES  6-1 


vii 

LOCKHEED  MISSILES  &  SPACE  COMPANY 


2-57-62-2 


Appendix 


A  TABLES 

B  SUPPLEMENTARY  NUCLEAR  WEAPON  EFFECTS  STUDY 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


I 

I 

I 

I 

I 

t 

I 

1 

I 

i 

I 

! 

1 

I 

I 

I 

I 

I 

I 


2-57-62-2 


ILLUSTRATIONS 


Fifure  Pagg 


4-1 

HuMd's  Case 

4-7 

4-2 

Hund's  Case  "b" 

4-7 

4-3 

-  X  Eg  Band  of  (Schumann-Runge  System) 

4-13 

4-4 

-3--  3  +  - 

B  II  -  A  2  y  Band  of  (1st  Positive  %Stem) 

4-15 

4-5 

3  3 

C  Ily  -  B  Ilg  Band  of  N2  (2nd  Positive  System) 

4-18 

4-6 

B  -  X  Band  Of  (1st  Negative  System) 

4-23 

4-7 

B  ^11  -  X  ^11  Band  of  NO  {0  System) 

4-25 

4-8 

A  ^E  -  X  ”2  Band  of  NO  (y  System) 

4-29 

5-1 

Effect  of  Haif-Width  Variation  on  Optical  Properties  of 

Og  for  V  *  31500  cm^^  and  T  ^  2000*  K 

5-4 

5-2 

Effect  of  Half-Width  Variation  on  Optical  Properties  of  Og  for 

Og  for  =  47500  cm'^  and  T  ^  2000*  K 

5-5 

5-3 

SACHA  Code  Average  Transmission  Result  for  Og  (0  5=  0.75) 

5-6 

5-4 

SACHA  Code  Average  Transmission  Results  for  Og  (cr  ~  1. 0) 

5-7 

5-5 

SACHA  Code  Average  Transmission  Results  for  Og  (or  =  2.0) 

5-8 

5-6 

Average  Transmission  Through  1  cm  of  0„  (calculated  from 
results  of  Meyerott,  Sokoloff,  and  Nicholll,  1959).  Electronic 
f-number  and  populations  normalized  to  agree  witii  SACHA  results 

5-10 

ix 

LOCKHEED  MISSILES  ft  SPACE  COMPANY 


Section  % 

INTRODUCTION 

This  report  surnwarizes  recent  progress  made  toward  the  ^eoretical  determination 
of  the  speciTai  absorption  coefficient  for  heated  air  in  the  temperature-density  region 
where  discrete  molecular  transitions  are  dominant.  Previous  investigations  at  this 
laboratory  have  been  made  along  similar  or  related  lines,  but  the  molecular  spectra 
were  not  considered  in  such  detail.  Meyerott,  Sokoloff,  and  Nicholls  (Ref.  l)  have 
calculated  the  Spectral  contribution  to  the  absorption  coefficient  of  air  for  the  opti¬ 
cally  thin  case,  while  Armstrong  et  al.  have  treated  the  similar  atomic  problem 
(Refi  2),  The  regions  in  the  temperature-density  plane  where  molecular  phenomena 
might  be  ejmccted  to  make  a  sizeable  conu^ibution  are  pointed  out  in  a  recent  article 
by  Armstrong,  Sokoloff «  Nicholls,  Holland,  and  Meyerott  (Ref.  3). 

In  the  present  study,  six  molecular  systems  will  be  considered:  the  Schumann-Rimge 
system  of  O2  ,  the  first  and  second  positive  systems  of  N2  >  the  first  negative 
system  of  N^  i  mid  the  beta  and  gamma  systems  of  NO  (nitric  oxide).  Future  calcu¬ 
lations  will  probably  include  any  continuum  as  well  as  0  photo-detachment.  Other 
effects  may  be  considered  as  the  tempertUure-density  regimi  is  extended  toward  that 
of  atomic  phenomena. 

After  necessary  preliminaries  are  developed  in  lection  2,  the  spectral  absorption 
coefficient  is  m^oduced  in  its  general  form  (for  molecular  transitions)  in  lection  3. 
lection  4,  after  a  review  of  spectroscopic  notation,  specializes  the  expressions  of 
lectimi  3  to  the  six  systems  of  interest  here.  Machine  calculatimis  and  results  are 
then  summarized  in  Sectimi  5,  in  which  a  preliminary  transport  code  is  also  discussed. 
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I  Section  2 

*  RADIATION  TRANSPORT 

®  2. 1  THE  EQUATION  OF  RADIATION  TR^SFER 


Each  volMme  element  in  a  mass  of  high  temperature  gas  is  Gimultaneously  emitting 
and  absorbing  radiation  of  all  frequencies.  It  is  the  balance  between  these  processes 
which  determines  the  local  radiation  flux  within  the  gas  as  well  as  the  flux  which 
escapes.  Let  the  flux  of  radiation  per  imit  frequency  mterval  about  u ,  per  unit  solid 
angle  in  the  direction  fl  and  per  unit  area  normal  to  the  direction  fl  ^  be  denoted  by 
I^(fl)  .  Similarly,  denote  the  energy  radiated  by  a  unit  mass  of  gas  per  unit  frequency 
per  lUiit  solid  angle  by  ji/  (fl):  ^  If  the  absorption  coefficient  of  the  gas  for  radiation  of 
frequency  u  is  k  cm^gm^^  ,  the  flow  of  radiation  throu^  the  gas  is  governed  by  the 
equation  of  radiation  transfer 

i 

p  <2.i) 

where  s  denotes  length  measured  in  the  direction  fl  . 


R  the  gas  is  in  local  fliermodynamic  equilibrium,  i.  e. ,  R  at  each  point  in  the  gas  a 
temperature  T  can  be  defined  such  that  the  characteristics  of  the  gas  m  the  neighbor¬ 
hood  of  the  point  are  those  of  a  gas  in  equilibriuni  at  that  ten[4>eraturei  the  emission 
coefficient  is  the  sum  of  a  spontaneous  and  induced  emission  term: 


'  IcT 

i  =  K  B  (T)  +  e  K  I  (2. 2) 

H  V  v'  '  V  V  \  •  / 

where 

I 
I 
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Md  is  tile  Piarick  distribution  fmction,  whieh  is  the  value  would  have  within 
ah  enclosure  in  thernaodynamie  equilibrium  at  temperaUire  T  .  The  following  dis- 
Gussion  is  limited  to  gases  in  local  thermodynamic  equilibrium. 


Substitutioh  of  Eq.  (2. 2)  into  Eq.  (2. 1)  yields 


dl 


(b  -  4- 

fl  as  p\v  vT 


-M- 


From  this  equation  it  is  clear  that  one  parameter  characteristic  of  the  gas  suffices  to 
specify  completely  its  radiative  properties,  namely 


t 


f 


PK. 


(2.5) 


the  absorption  coefficient  of  the  gas. 

In  many  problems  the  spectral  distributiQn  of  the  radiation  is  not  of  primary  concerni 
For  these  cases  an  appropriately  defined  mean  value  of  the  absorption  coefficient  is 
useful^  The  manner  in  which  the  mean  value  is  calculated  depends  on  the  character¬ 
istics  of  the  problem  under  consideration.  A  limiting  case  of  interest  is  that  of  mi 
optically  thin  sample  of  gas,  i.e, ,  a  sample  whose  dimensions  are  small  compared 
with  the  mean  free  path  of  radiation  in  the  gas.  Consideration  of  this  case  leads  to 
the  definition  of  the  Planck  mean  absorption  coefficient.  At  the  opposite  emxeme  is 
the  case  of  an  optically  thick  gas  sample,  which  is  conveniently  analyzed  in  terms  of 
the  Bosseland  mean.  These  two  limiting  cases  are  discussed  below.and  the  respec¬ 
tive  mean  absorption  coefficients  are  derived. 

2. 2  THE  PLANCK  AND  ROSSELAND  MEAN  ABSORPTION  COEFFICIENTS 

Consider  the  radiation  from  an  isolated  thin  slab  of  gas.  The  total  mmiocluroniatic 
radiaticii  flux  per  unit  area  leaving  one  face  of  the  slab  is 

=  yi^cosOdn  (2.6) 

2ff 
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where  ©  is  measiured  from  the  normal  to  the  slab.  Multiplying  Eq.  (2. 4)  on  both 
sides  by  dt2  and  integrating  over  the  hemisphere  of  outward  directions  yields  the 
equation 


cos  6  dQ 


2V  /i  B 

V  V 


(2.7) 


where  x  is  dist^ce  normal  to  the  slab  face.  From  Eq.  (2. 6)  and  the  definition  of  the 
mean  outward  intensity  (Ref.  4,  Equation  3.43) 

=  A 

2Jr 


Eq.  (2. 7)  may  be  rewritten 


(2.8) 

Our  thin  slab  approximation  now  consists  in  neglecting  relative  to  B  ,  and  taking 
cl^^ '  .j. 

as  the  ratio  of  finite  increments  of  emitted  flux  5F  and  slab  thickness  5x  . 

With  these  approximations, we  obtain  from  Eq.  (2. 8)  for  the  integrated  flux  or  total 
radiant  energy  endtted 

Sf  =  27r  fix  /  P,  B  di/  (2. 9) 

V  V 

The  corresponding  radiation  emitted  by  a  perfect  radiator  in  equilibrimn  at  tempera*^ 

4 

ture  T  is  oT  ,  The  emissivity  c  of  &e  gas  is  tiOe  ratio  of  its  total  radiatimi  to 
that  emitted  by  a  black  bocty  at  the  same  temperatuure;  hence 

2ir  fix  /#  B  dv 

J  -V 

e  ^  (2.10) 
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The  quantity  ^  ,  the  emisslvlly  per  unit  length,  is  usually  quoted.  It  can  be  written 


-L. 

stV  " 


B^(T)d>.  •  2S^(T) 


(2. 11) 


where 


WTB- 


(T) 


(2. 12) 


is  the  Planck  mean  absorption  coefficient,  sometimes  referred  to  as  the  emissimi 
mean  absorption  coefficient. 


Now  consider  the  case  where  the  dimensions  Of  the  gas  sample  are  large  compared 
With  the  mean  free  path,  so  that  radiation  cannot  move  freely  from  one  part  of  tiie  gas 
to  another.  To  first  sq>proximati(m  we  take 


(2.13) 


Substituting  for  in  the  transport  equation,  in  the  derivative  only,  gives 


I,  =  B  -  s  ^ 
<*S  lA^  ds 


(2.14) 


and  re-sid>stitution  wotdd  produce  higher  order  terms.  Shice  is  isotropic,  the 


e^qpressitm  for  the  x  component  of  the  flux  reduces  to: 


■//■ 


Fjj  -  I  I  Ij,  cos  9  dn  dv 


COS' 9  sin  9  d9  d0  dv 


.  _  4£  J_  ^  f  ^ 

3  jS^  dx  J  dT 


,  4r  1  dBf 
’lx 


(2.W) 
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where  ©  is  measured  from  the  x  axis,  B  (T)  is  the  integrated  Planck  intensity, 
and 


X 


(2. 16) 


defines  ,  the  Rosseland  mean  absorptidn  coefficient.  Related  quantities  are  the 

rt 

Rosseland  mean  opacity,  defined  as  tiu  =  ^  iSnd  the  Rosseland  mean  free  path 

A  ^  1  ^  ^ 

» '  h  • 


Equation  (2. 15)  may  then  be  written  in  the  form 

f  ^  ^  Ag  Hit) 


(2, 17) 


The  condition  for  validity  Of  Eq.  (2. 15)  and  (2*  17)  is  that  the  fractional  variation  m 


temperature  be  small  in  a  distance  of  one  mean  free  path  for  all  frequencies  of 


interest;  hence 


(2. 18) 


2. 3  IMPORTANCE  OF  LINE  WIDTHS  IN  RADlA'nON  TRANSPORT 

The  in^ortance  of  Ime  effects  in  radiation  transport  can  be  illustrated  most  easily  in 
the  two  <^posite  limits  of  die  optical  depth  for  which  simple  solutions  to  the  transport 
equation  exist,  namely,  in  the  limit  of  mi  optically  diin  sample  and  in  the  limit  of  an 
optically  thick  sample.  These  two  cases  were  discussed  m  Section  2. 2.  Consider 
first  the  optically  dense  case  where  the  Rosseland  mean  absorption  coefficient  []Eq. 

(2. 16)3  becomes  a  useful  parameter.  The  Rosseland  mean  is  an  inverse  mean  and 
thus  depends  critiically  on  the  minimum  values  M ( v)  ,  the  so^alled  "windows,  " 
Conversely,  it  is  relatively  insensitive  to  the  intensities  at  the  peaks  of  the  lines. 
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especi^ly  for  strong  lines.  The  most  importaiit  characteristic  of  the  strong  line  is, 
in  the  optic^ly  thick  case,  the  width  over  which  it  effectively  blacks  out  the  spectrum. 

In  the  case  of  an  optically  thin  sample  the  useful  parameter  is  the  Flmick  mean  absorp¬ 
tion  coefficient  [Eq.  (2. 12)*],  Assuming  the  Imes  to  be  sufficiently  narrow  so  that 
the  factor 

hi' 


does  not  change  appreciably  across  the  width  of  a  line,  we  see  that  the  line  shape  is 
clearly  unimportant.  Each  line  can  be  considered  a  delta  function  and  the  frequency 
integration  performed  immediately.  It  follows  that  the  Planck  mean  is  an  additive 
mean  and,  as  such,  can  be  used  to  estimate  the  gross,  overall  contribution  of  the  lines 
to  the  total  absorption  coefficient. 


The  evaluation  for  the  intermediate  case  cmmot  be  simply  determined,  but  one  must 
resort  to  a  detailed  transport  calculation  to  determine  the  sensitivity  to  the  widths  and 
sh^es  of  the  lines.  A  calcvdation  has  been  carried  out  to  estimate  the  effect  of  line 
widths  on  a  transport  problem.  This  calculation  will  be  described  in  Section  5. 
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Section  3 

THE  LINE  ABSORPTION  COEFFICIENT  OF  A® 

3. 1  EXPRESSION  FOR  THE  ABSORPTION  COEFFICIENT  OF  A  SINGLE  LINE 


We  first  derive  an  ej^ression  for  the  absorption  GOefficient  for  a  single  ''pure"  linej 
that  is,  a  line  arising  from  a  tr^sition  between  non-degenerate  electronic  levels  or 
sub-levels.  *  In  the  following  it  will  be  assumed  that  the  reader  has  a  reasonable 
familiarity  with  conventional  speGtrosc^ic  notation  and  widi  the  theory  of  the  spectra 
of  diatomic  molecules  on  a  level,  say,  with  that  in  Herzberg's  text.  Spectra  <rf  Diatomic 
Molecules  (Ref.  5).  For  the  reader's  ccmvenlence,  however,  we  have  mcluded  some 
background  material  in  Section  4. 

The  line  absorption  coefficient  for  a  transition  from  a  lower  state  (n" ,  v" ,  J") 

to  an  iQ>per  state  ( n' ,  V  *  J' )  is  given  by 


where 


M(v)  *  N^„^„J„  ^*'n"n',V"vSJ"J' 

n  t  ,1  T"  "  Bohr  freqnenCy  (measured  in  cm"'^)  of  the  line 
^n"v"J"  ^  particle  density  in  the  lower  state  (n"  ,  v"  ,  J") 
B  „  ,  „  ,  T„  T,  ’  Einstein  coefficient  for  induced  absorption 

11  llfV  V>v  V 

F  ( )  =  a  line  shape  factor  such  ttiat  /  F  {u)dv  =  i  . 


(3.1) 


♦By  svdt)=level  we  simply  mean  one  of  the  components  of  a  degenerate  electronic  state, 
for  example,  one  of  the  25+  i  components  of  a  S  state  with  spin  s  or  one  of  the 
conmonents  of  a  A-type  doublet.  By  electronic  degeneracy  we  mean  ttie  spin  and 
orbital  degeneracy  that  would  exist  m  the  absence  of  nuclear  motion  (rotation  mid 
vibration),  spin-orbit  interactions,  and  electron  spin  interactions,  (We  ignore  com¬ 
pletely  any  interaction  of  the  nuclear  spins  with  the  rest  of  die  molecule. )  The  2J  + 1 
spatial  degeneracy  associated  with  raolecidar  rotation,  which  rigorously  exists  in  an 
isotropic  envirmiment,  is  never  counted  in  figuring  electronic  degeneracy. 

♦♦Doiflile  and  sinsde  primes  will  be  used  to  denote  lower  and  igiper  states,  respectively. 
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Here  n  is  used  to  denote  a  non-degenerate  electrcffiie  level  or  sub-level,  whUe  v 
and  J  denote  vibrational  and  rotational^  levels,  respeotively. 


The  Einstein  coefficient  j„j,  is  given  by 

E 


4^ - — 


M"M' 


'WW 


®n''n',v"v',J'’J'  ~  ».2 
3n  c 


2J"+  1 


where 


rM"M'  ^  /•  u  m'  df 

J  n"v"  J"M''  “^n'  v'J'M' 


(3.3) 


is  the  matrix  element  of  the  electric  moment  operator  R  =  H  +  H-  of  die  electrons 

e  n 

and  nuclei.  M"  and  M'  are  azimuthM  quantum  numbers  numbering  the  (spatially) 
degenerate  rotational  levels  of  the  lower  and  upper  states,  respectively.  The  summa¬ 
tion  in  Eq.  (3. 2)  is  over  all  possible  combinations  of  the  rotational  sub-levels  of  the 
lower  with  those  of  die  upper  state.  If  we  assume  the  Born-Qppenheimer  approKimation 
to  be  valid,  the  molecular  wave  function  may  be  separated  into  electronic,  vibrational, 
and  rotational  parts 

^  (3.4) 

The  vtorational  eigenfimction  depends  on  the  internuclear  distance  oidy.  The 
electronic  eigeiffunction  4'  depends  on  the  electron  coordmate  and  slightly  on  the 
internuclear  distance  as  a  parameter.  The  rotational  eigenfimction  4^^  depends 
on  the  nuclear  coordinates  and,  in  the  general  case,  also  on  the  electron  coordi” 
nates ,  inasmuch  as  the  quantum  number  J'  here  refers  to  the  total  angular  momentum 
resulting  from  nuclear  rotation  and,  if  present,  electron  spin  and  orbital  angular 
momentum. 
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By  substituting  Eq^  (3.4)  Into  Eq.  (3  3).  it  is  readily  Shown  that 


9  .  .9 


where 


t  ^  ti  I 


q(v’v")  ^  I  %  dty 


Here,  Sj„  is  a  luie  Strength  factor  and  is  that  part  of  E  j;  R  |  that  depends  on 
J"  (and  also  on  A  ,  the  component  of  the  electronic  orbits  angular  momenhim  along 
the  internuclear  axis,  whenever  A  ^  O).  Sj  is  usually  called  ^e  Honl-Londen 
intensity  factor.  HonULondon  factors  for  all  the  important  types  of  electronic  transi¬ 
tions  are  availt^le  In  the  literature.  is  the  electronic  transition  moment  which, in 
general, depends  on  the  internuclear  separatipn,  q(v'v")  is  the  so-called  Franek- 
Oondon  factor  for  the  (v ,  v")  band. 


Equation  (3. 2)  then  becomes 


n''n',v"v',J"J' 


I  «el  prn 


The  particle  density  (n”,  v”,J")  is  given  by 


N  =  N 

%"v"J"  total 


(2J”  +  1)  Wj  exp  («E^„^„j„At) 


where  is  the  total  particle  density  for  all  levels,  ®nTryi»  j??  is  the  energy  ^  tiie 

level,  Uj  is  the  nuclear  spin  statistical  wei^t  for  this  level,  and  is  the  total 

partition  funcHon  for  the  system. 
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The  total  partition  function  may  be  written  as  a  sum  of  contributions  from  each 
electronic  state: 


n 

where  n  denotes  the  electronic  state.  Taking  into  account  elecironic  (orbital  +  spinl, 
vibrationalj  rotational,  and  nuclear  (spin  only)  degrees  of  freedomi  we  may  write 
in  the  form 


“  ^electronic  ^vib-rot. 


Q, 


nuclear 


where 

hcv 

‘’electronic  “  “A<“  ♦  «  >  “A  =  1  “  A  c  0 

=  2  a  A  ¥  0 


00 


O 


(2Ij^+l)(2J^+ 1) 

^nuclear  *  c 


(3. 11) 


(3. 12) 


(3.13) 


(3. 14) 


(3. 15) 


In  these  formidas,  v  is  the  energy  of  the  lowest  vibrational  level  of  electronic  state 
oo 

n  above  that  of  the  grotmd  state,  Gq(v)  is  the  vibrational  term  referred  to  the  lowest 
vibrational  level  as  zero,  Fy(J)  is  the  rotational  term  for  the  vibrational  level 
and  is  the  correspondhig  rotational  constant,  is  the  statistical  wei^t  for 
orbital  angniar  momentum  (A)  about  the  intemudear  aids,  S  is  the  totsd  ^in  of 
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the  eleclrdnSj  aiid  are  the  nuclear  spins  of  ^e  two  fluelei  a  and  b,  respec¬ 
tively,  and  cr  is  a  symmetry  number  havmg  the  value  2  for  homonuclear  molecules 
and  1  for  heteronuclear  systems. 


An  approximate  formula  for  Qyjjj-  jiQt  fs  somewhat  easier  to  use  than  Eqs.  (3.13) 
and  (3. 14)  and  yet  entirely  adequate  for  our  purposes  has  been  given  by  Bethe  (Ref.  6) 


erlvmef.  7i- 


*vib-rot. 


1  -  ejm  (“  1.4388  oi 


1.4388 


(1+YT) 


(3. 16) 


where  y  ,  Bethe 's  correction  factor  for  anharmonicity  and  non-rig^dily,  has  the  value 


,  /2(o  X  O;^  8B  \ 

^  - 1 - 

1.4388  w.  \  0)  B  w  I 

0  \  Q  o  oy 


(3. 17) 


Here  0  and  B  are  the  vibrational  and  rotational  constants  of  Ihe  v  =  o  vibrational 
0  0  -  - 

level,  0  X.  is  the  first  anharmoniq  constant  and  a  is  the  interaction  constant  Iw 
0  0  .  o 

cov^ling  between  rotation  and  vibration  (Ref.  5). 

Rehurning  now  to  Eq.  (3. 9)  we  rnttitiply  numerator  and  denominator  by  as  given 
by  Eqs.  (3. 11  ^  3. 15)  and  obtain  after  some  manipulatimi 


XT  Sa"  I 

^n"v"J''  *^total  (3S+ 


mclear 


(3. 18) 


(2J"+1)  expj^-  j  G^(v")  +  F^„(J")(  hc/kTj 
^vib-rqt. 


^ere  use  has  been  made  of  the  relation 


i"v"J"  ^  [‘'oQ 


(3. 19) 
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Si^stitutiflg  Eq,  (3.  8)  for  j.i  jr  (3. 18)  for  N^Mynji. 

Eq.  (3. 1)  we  obtain  the  following  rather  complicated  e}q)ression  for  the  line  absorption 


coefficient: 


M(v) 


^  [^total  ^n"] 


8Jf^  u>, 


V%n-.v- v'. J". ^ 
exp  ]g^(V^)  +  Fy,r(J*’)} 


q(v*.v")Sj„ 


Q, 


F(v) 


vib-rot. 


(3.  20) 


^.th 


where  F^i>  ^  Qn-Z^total  is  the  fractional  population  of  the  n"  electronic  state. 
By  introduction  of  the  notation 


8#^  L, 


Vn».v"v*,J"J' "  3hc°  w^l2S+l)  Qnuclear  v" V, J»J' l^e I  )®j'‘ 

(3.21) 


=  [Gq{v’')  +  f”„(J")]  hcA 


(3. 22) 


where  L  ???  2.6875  *  10^^  particles/cm^  is  Loschmidt's  Numbert  Eq.  (3. 22) 
0 

becomes 


p(v)  ^ 


N  P 
^total  -  n" 

^o  ^vtt)-rot. 


Note  that  the  quantities  H  and  E  are  characteristic  of  die  isolated  moleeide,  that  is, 
do  not  depend  on  the  temperature  and  densi^.  Note  that  the  bracketed  term  m  Eq. 

(3. 23)  is  dimensionless  and  that  H  has  die  dimensicais  cm  .  Since  the  sh^e 
factor  F(i/)  has  dimensiims  cm  (see  below),  I4(i')  has  the  dimension  cm  as 
required. 
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Although  Eq*  (3^  23)  femalns  valid  for  arbitrary  form  factor  F(v  )  ^  in  all  Of  Our 
work  we  have  assumed  a  Lorentz  line  profile,  viz. , 

P(.)  .  i  .  ^  (4.24, 

O  +  yHyl^  Jll  Jty 

where  cr  is  the  line  half -width  at  hMf-maxlmum. 

3.  2  OVERLAPPING  AND  COMPOSTE  LINES 

The  absorption  coefficient  ti(v)  at  v  due  to  severai  overlapping  lines  is,  of  course, 
the  sum  of  the  separate  line  absorption  coefficients: 

M(v)=£u.(v)  (3.25) 

i  ^ 

where  i  labels  the  lines  and  the  sum  runs  over  all  the  lines  whose  profiles  overlap 
signilicantly  at  v  . 

Fine  Structure  may  frequently  be  ignored  without  compromising  the  accuracy  of  the 
calculation.  In  these  cases  the  intensity  factor  of  die  unresolved  line  may  be  obtained 
by  combining  the  intensity  factors  of  the  imresolved  components.  An  example  of  the 
consfructiott  of  sueh  an  intensity  factor  will  be  included  in  Section  4. 5. 


I 
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Section  4 

THE  LINE  FREQUENCIES  AND  LINE  INTENSITIES  FOR  ADR 


4. 1  SPECTROSCOPIC  QUANTUM  NUMBERS 


In  a  diatoniic  molecule^  the  synametry  of  the  field  in  which  the  electfons  Move  is  axial 
about  the  line  connecting  the  two  nuclei,  with  the  coMponeat  of  ^e  orbital  angular 
MOMentum  of  die  electrons  about  the  interauclear  axis  being  a  constant  of  the  motion, 
The  orbital  angular  momentum,  L  ,  precesses  about  this  axis  with  constant  eompo= 
nent  Mj^ii  ,  where  is  restricted  to  the  values  Mj^  =  L,  L-1,  L“2,  ‘  ,  “L  . 

Since,  in  diatomic  molecules,  states  differing  only  in  the  sign  of  have  the  same 
energy  (are  degenerate),  it  is  appropriate  to  introduce  the  quantum  nunober 


A.  =  Mj^l  ;  s  0,  1,  2, 


If  L  precesses  very  rapidly,  its  meaning  as  angular  momentum  noay  be  lost,  but 
and  A  remain  well  defined: 

When  A  =  0,  1,  2,  3,  •  •  •  ,  the  corresponding  molecular  state  is  designated  a 

S  ,  n  >  A  ,  •  state,  analogous  to  the  mode  of  designation  for  atoms.  By  con» 

vention,  Greek  letters  refer  to  the  components  of  electronic  angular  momenta,  wMle 

the  corresponding  lower  case  letters  used  for  atoms  refer  to  the  electronic  angular  momenta 

themselves.  11 ,  A  ,  ^  ,  "  *  states  are  doubly  degenerate  since  can  have  tiie 

two  values  +  A  and  =  A  .  S  states  are  non-^degenerate. 

Just  as  for  atoms,  the  spins  of  the  electrons  form  a  resultant  S  ,  the  corre^onding 
quantum  nunyher  s  being  integral  or  half  integral  as  determined  by  whether  the  total 
number  of  electrons  in  the  molecule  is  even  or  odd.  The  component  of  s  along  the 
axis  of  symmetry  is  designated  s  (analogous  to  atoms) ,  and  this  must  not  be  confused 
with  the  Z  which  means  A  ^  0  ,  The  values  allowed  Z  are  S,  S-l,  S-2,  •  •  •  ^S, 
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ii  ,  2S  4-  1  deferent  values  are  possible.  The  quantum  number  S  can  be  either 
positive  or  negative,  except  it  is  not  defmed  when  A  =  0  (for  S -states). 


The  total  electronic  angular  momentum  about  the  internuclear  axis  is  obtained  by  an 
algebraic  addition  of  a  and  S.  The  resultant  quantum  number  is  given  by  -I  A  +s|. 
tt  A  ^  0  ,  there  are  2S  +  1  values  of  A  +  2  for  a  given  value  of  A  .  The  inter¬ 
action  of  S  with  the  ma^etic  field  produced  by  A  gives  rise  to  a  multiplet  sb?ueture. 
"The  term  symbor  carries  2S^+  iras^  a  t^superscrlpr  and~A  +  s  as  a  ri^t  subscript, 
for  example  the  triplet  a  would  have  components  a  ^  ,  A  g  ,  and  a  ^  ,  Even  for 
2  states  (there  is  no  magnetic  field  m  the  direction  of  the  internuclear  axis)  the  multi¬ 
plicity  is  said  to  be  2S  +  l  even  though  no  actual  splitting  may  occur  for  non-rotating 
molecules. 


4. 2  SYMMETRY  PROPERTIES  OF  THE  EIGENFUNCTIONS 


Under  the  coordinate  transformation  X|^— *«3^,  yj^-=^-yj. ,  ,  the  wave 

function  for  any  atomic  system  either  ch?mges  sign  or  remains  unchanged  in  sip. 

For  a  rigid  rotator,  a  reflection  at  the  origin  is  accomplished  by  the  transformation 
0  ir-O,  IT  +  ^  ,  It  turns  put  Hiat  the  rotator  wave  function  remains  unchanged 
for  even  J  and  changes  sip  for  odd  J  ^  In  the  case  of  the  symmetric  top  (A  **  0> , 


each  rotationeil  level  is  doubly  degenerate.  Linear  combinaticms  that  change  sip  or 
remain  unchanged  under  reflection  can  always  be  chosen,  however, 


The  rotational  levels  of  a  diatomic  molecide  are  classified  according  to  their  behavior 
voider  reflection  of  the  total  wave  function  (not  just  the  rotational  wave  function  alone) . 
A  rotational  level  is  called  positive  or  neptive  depending  on  whether  tiie  total  wave 
ftmction  remains  unchanged  or  changes  sip  for  reflection  at  the  origin.  To  a  first 
approidmatimi,  the  total  wave  function  for  a  diatomic  molecule  can  be  written 

♦  (4.1) 

6  V  r  '  f 


Since  it  depends  only  on  the  magnitude  of  the  internuclear  distance,  remains 

V 

unchanged  under  reflection.  Consider  now  two  cases,  A  ^  0  and  A  ^  0  . 

4-2 
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I  4.  2. 1  Case:  A  =  0 

■  ^  %  remains  unchaaged  under  refleetion,  the  pmfity  depends  only  on  Md  the 

*  rotational  levels  are  +  or  -  ,  depending  on  whether  J  is  even  or  odd.  (See  below.  )i 

■  ^  %  changes  sig^i  under  reflection,  however,  the  total  wave  function  remams 


unchanged  for  odd  J  and  changes  sign  for  even  J  .  In  case  odd  J  levels  are 


Energy  Level  Diagram 
(I  is  the  operator  which  replaces  r  by  -r  ) 

4. 2.  2  Case;  A  ^  0 


In  diis  ease  for  each  value  of  J  there  are  positive  and  negative  rotation^  levels  of 
equal  energy.  (See  above. )  The  levels  are  drawn  as  if  they  were  degenerate.  Actually,, 
when  account  is  taken  of  the  interaction  between  elecn?onic  motion  and  rotation,  there 
is  a  slight  splitting  of  the  degeneracy.  Note  that  the  J  ^  Q  level  cannot  occur  for 
A  -  1,  since  J  here  is  the  total  (rotational  plus  orbital)  anpilar  momentuni. 

It  must  be  emphasized  that  the  positive  or  negative  character  of  the  rotational  levels 
depends  (mly  on  the  rotational  symmetry  of  the  system,  i,  e. ,  equivalence  of  r^t  and 

4*3 
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left-handed  coordinate  systems,  £md  has  nothing  to  do  with  the  approximate  represen¬ 
tation  of  the  total  wave  funetion  given  by  Eq.  (4. 1).  This  approximation  was  introdueed 
only  for  the  purpose  of  determining  which  levels  are  positive  and  which  are  negative, 
it  will  be  shown  below  that  positive  levels  combine  only  with  negative  levels  under 
dipole  radiation  and  vice  ver  sa.  Symbolically 

.  I  |-  I  .  f- - - 

Clearly^  this  is  compatible  with  the  selection  rule  A  J  -  ±  1  for  the  simple  rotator  or 
A  J  -  ±1  for  the  symmetric  top  rotator  with  A  ^  0  . 

The  <mposlte  selection  rule  holds  for  the  Raman  effect;  +  terms  combine  only  with  + 
and  -  Only  with  -  ,  or  symbolically 

which  is  consistent  with  the  selection  rule  A  J  ^  6 ,  ±  2  for  the  simple  rotator  or 
Aj  =  0,  ±1,  ±2  for  a  symmetric  t(^  with  A  0  . 


Selection  rules  for  dipole  radiation  depend  on  integrals  of  the  form 

M  dr 
y  n  m 

Now  M  -  =1M  .  If  both  'if  and  are  positive  or  both  negative  under  reflection, 
the  integral  will  change  sign.  This  is  only  possible  if  the  integral  is  identically  zero, 
i.e.,  such  transitions  are  forbidden.  On  the  other  hand,  transitions  between  positive 
and  negative  levels  are  allowed  since  the  sign  of  the  integral  does  not  change  under 
reflection. 

The  Raman  rule  follows  in  the  same  way  by  noting  that  the  polarizability  is  unchanged 
under  reflection. 
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For  homonuclear  molecules  rotational  levels  may  be  further  Glassified  according  to 
their  behavior  mder  interchaage  of  the  two  identical  nuclei.  In  the  same  way  as  one 
shows  thatjfor  the  exchange  of  two  electrons,  the  total  wave  function  must  either 
change  sipi  or  remain  unGhanged  in  sign,  we  can  show  that  for  an  exchange  of  nuclei 
the  total  wave  function  either  remains  unchanged  in  sigh  or  changes  only  its  sighi  in 
which  case  the  state  under  consideration  is  said  to  be  symmetric  or  antisymmetric  in 
the  nuclei,  re^ectively.  We  can  also  show  that  in  a  given  eleetronie  state  either 
the  positive  rotational  levels  are  symmetric  and  the  negative  are  antisymmetric 
throughout^  or  the  positive  rotational  levels  are  antisymmetric  and  the  negative  are 
symmetric  throughout^  As  an  example  consider  the  case  A  =  0 .  Clearly 
there  are  four  cases: 


J 

4.-.r:--r  -  +s 

S'*" 

g 

3  r;T7r---  - 

2  -  -  -  - -  +8 

1--  -  -a 
0 - - - — +s 


J 

4  ■  +a 

t' 

u 

3  -  - — .-S 

■  »a 

1  ■  -■ — s 

0— -^+a 


J 

u 

3— — ^+a 
2 —  -S 

1--  ■  ■-*a. 

0— s 


J 


4  — — — — ^-a 

s 

2-— - - — a 

1 . . —  *s 


♦ 

e 


♦ 


4-  = 

e 


A  -  0 

K  the  nuclei  have  zero  nuclear  spin  (or  if  the  interaction  between  nuclemr  spin  with  &e 
rest  of  Uie  molecule  is  neglected),  there  is  a  rigorous  prohibitimi  of  intercombinatimi 
between  symmetric  and  antisymmetric  states: 
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This  holds  for  transitions  brought  about  in  any  ways  e.  g. ,  collisions,  Eaman  effect, 
dipole  radiation,  etc.  The  proof  of  this  follows  trivially  from  the  form  of  the  operator 
for  the  type  of  transition  involved.  Physically,  operators  for  dipole  radiation ^  colli¬ 
sions^  Raman  effect,  etc,  ^  must  be  symmetric  in  the  two  nuclei.  Since  transition 
probabilities  are  governed  by  integrals  of  the  form 


it  follows  the  and  ^  must  transform  in  the  same  way  for  interchange  of  die 
n  m  ^  a 

nuclei. 

An  immediate  consequenee  of  the  non-intereombmation  rule  is  that  there  are  no 
vibration -rotation  or  rotational  spectra  for  homonuclear  diatomic  molecules  since  two 
rotational  levels  for  which  the  selectimi  rule  A  J  =  ±  1  is  fulfilled  have  opposite 
symmetry  in  the  nuclei.  Raman  effect  transitions  can  occur,  however,  with  A  j  =  0, 
±2. 

For  quadrupole  radiation,  AJ  -  0,  ±1,  ±2  .  It  is  seen  that  AJ  -  0,  ±2  is  compat¬ 
ible  with  aw  I  ‘8  and  therefore  this  type  of  radiation  can  occur  in  homonuclear 
molecules. 

4.3  COUPUNG  OF  THE  ANGULAR  MOMENTA 

There  are  mutual  interactions  between  the  electronic,  vibrational,  and  rotational 
motions  of  a  diatomic  molecule.  The  electronic-vibrational  interaction  is  taken  into 
account  by  choosing  the  vflbrational  levels  to  fit  the  electronic  potential  curve;  the 
vibrational -rotational  mteraction  is  accounted  for  by  considering  the  molecule  to  be  a 
vibrating  rotator.  Mutual  mteraction  of  die  electronic  and  rotational  motions  must  be 
cmisidered  also;  this  interaction  depends  on  which  quantum  numbers  distinguish  the 
rotational  levels  in  different  types  of  electronic  states,  and  on  the  energy  dependence 
on  the  quantum  numbers, 
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Angulo  momenta  in  the  molecule  combine  to  form  a  resultant  desi^ated  £  .  The 
spin  S  ,  nuclear  rotation  N  ,  and  orbital  angular  mementum  A  are  coupled  in 
dMferent  ways  to  form  •  3  The  modes  of  coupling  were  classified  by  Hund  and  are 
called  Hund's  case  "a,"  "b,"  "c/’ —  etc.  Of  these  we  need  discuss  only  two, 
cases  "a"  and  "b." 

4. 3. 1  Hund's  Case  "a" 


In  Hund's  case  "a"  it  is  assumed  that  the  interaction  of  me  nuclear  rotation  with  the 
elecironic  motion  (orbital  *  spin)  is  very  weak,  whereas  the  electronic  motion  is 
strongly  coupled  to  the  internuclear  axis  (see  Fig.  4-1). 


Fig.  4-2 
Hund's  Case  "b" 


The  total  electronic  angular  momentum  il  is  then  well  defined,  n  and  N  coir^ine 
to  form  J  .  The  precession  of  L  and  S  about  the  mternuclear  axis  is  assumed  to 
be  very  much  faster  than  the  nutation  of  the  figure  axis,  il  is  integral  or  half-integral 
depending  on  whether  me  number  of  electrons  is  even  or  odd.  Since  J  cannot  be 
smaller  than  its  component  n  , 

J=0,  0+1,  0+2,  ••• 


Levels  with  J  <  0  do  not  occur. 
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4. 3.  2  Hand's  Case  "b" 

When  the  spin  is  not  coupled  to  the  ittternuclear  axis  at  all,  Q  is  not  defined  and 
Hund's  case  "a"  cannot  apply.  This  is  true  where  A  -  O  and  S  ^  0  and  May  be 
nearly  true  in  other  instances,  in  Mund's  case  "b,"  therefore,  the  spin  is  decoupled 
IcPM  the  internuelear  axis  or  is  very  weakly  coupled  to  it.  A  and  N  form  a  resultant 
K,  which  then  combines  With  s  to  form  a  resultant  J  (Fig.  4-2).  H  A  s  o  ,  it  is 
identical  with  N  and  is  therefore  perpendicular  to  the  internuelear  axis.  In  general 

K  may  have  the  integral  values  KsA,A+1,  A+2, .  .  The  possible 

values  of  J  are  (K+S)  ,  (K+s-i)  ,  (K+s-2)  ,  •••  [K-sj  ;  for  a  given  K,  Except 
when  K  <  S,  then,  each  level  of  a  given  K  consists  of  2S+1  components.  Again 
J  is  half-integral  or  integral  accOrdmg  to  whether  the  nundjer  of  electrons  is  odd  or 
even,  respectively. 

A  Type  Poubling 


In  Hund's  cases  "a"  and  "b>"  the  interaction  between  N  and  L  has  been  neglected. 
For  larger  speeds  of  rotation  this  interaction  rnust  be  taken  into  account  and  is  found 
to  produce  a  splitting  into  two  components  for  each  J  value  in  states  with  A  ¥  0 
which  are  dotmiy  degenerate  When  there  is  no  rotation.  This  splitting,  which  is  pre* 
sent  m  all  states  where  A  ¥  0  and  increases  with  increasing  rotation  ( J)  ,  is  called 
A -type  doublii^. 

Spin  Uncoupling 


It  may  ocem’  (and  frequently  does)  that  a  molecule  coupled  according  to  Hund's  case 
"a"  for  small  J  may,  with  increasmg  rotation,  go  over  to  case  "b,"  This  happens 
when  the  rotational  velocity  of  the  molectde  becomes  large  compared  with  the  pre- 
cessional  velocity  of  s  about  A  .  S  then  is  essentially  decoupled  so  diat  is  not 
defined  but  K  is.  This  process  is  called  spin  \mcovg>ling.  It  is  often  convenient  to 
formally  extend  the  quantum  number  K  back  to  small  rotation  values  even  thou^  it 
is  not,  strictly  speaking,  defined. 
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4.4  SELECTION  RULES  AND  IRANCHES 

J  ,  the  qMaiituffla  niutiber  of  the  total  angular  mOManttB,  Obeys  the  seleGtion  rule 
AJ  =  0  j  +1  j  -1  with  the  restrietioa  that  Aj  =  0  is  forbidden  when  A  =  0  in  both 
the  upper  and  the  lower  eieetronic  states.  This  seleetion  rule  forms  a  basis  for  elas* 
sifying  the  speotral  lines  into  distinct  series  called  branches,  J'  -  J”  +  1  eorres=‘ 
ponding  to  the  R^^branch,  J'  =  J''  -  1  to  the  P-branchi  and  J*  “  J''  to  tho  Q^brafich. 
Thus  the  rotational  contribution  to  the  frequency  of  a  transition  by  the  branches  can 
be  expressed  as 


P(J)  =  F'(J  “ 

1)  *  F"(J) 

Q(J)  ^  F'(J)  - 

F"(J) 

R(J)  =  F*(J  + 

1)  *  F"(J) 

The  above  are  called  main  branches  to  distinguish  them  from  ^ellite  branches .  For 

illustration,  assume  that  both  the  upper  and  the  lower  state  take  case  ''b"  coupling. 

Then  the  relevant  quantum  ntunbers  are  K  ,  S  ,  and  J  ,  The  transitions  with 

Aj  =  AK  are  the  main  branches,  and  those  with  AJ  *  AK  are  the  satellite  branches. 

The  satellite  branches  always  lie  close  to  the  main  branches  and  are  usually  much 
■  2  2  + 
weaker.  If,  for  example,  the  transition  takes  place  between  a  n(b)  and  a  2  state, 

with  J’  ^  j"  and  K'  K"  +  1  ,  then  this  transition  is  called  an  R<-£orm  Q  satellite 

branch  and  is  abbreviated  .  K  is  a  Q  branch  because  J'  ^  J”  and  r-form 

because  K’  =  K"  +  i  . 

Various  other  selection  rules  arise  from  the  symmetries  in  a  dynamical  system.  For 
example,  the  requirement  that  the  total  wave  function  for  a  molecule  cmi,at  most,change 
sign  for  an  inversion  of  the  coordinates  of  all  the  particles  leads  to  the  nde  that  in  the 
dipole  approximation  transitions  can  take  place  oidy  between  rotational  levels  of 
opposite  sign.  R  the  nuclei  of  a  diatomic  molecule  are  identical,  an  exchange  of  the 
nuclei  should  either  leave  the  total  wave  function  unchanged  or  change  its  sign  for 
symmetric  or  antisymmetric  states.  The  corresponding  gener^  selection  rule  pro* 
Mbits  transitions  between  syminetric  and  antisymmetric  states.  When  the  two  nuMei 
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are  not  identical  but  do  have  the  same  charge,  the  electronic  wave  function  should 
suffer  at  most  a  change  of  Sign  when  the  coordinates  of  the  electrons  are  changed  in 
sign.  If  the  electronic  wave  function  does  not  change  sipi,  the  total  molecitlar  wave 
function  is  "even"  or  "gerade."  Otherwise  it  is  "odd"  or  "ungerade."  Even  states  do 
not  eombine  with  even  states  i  nor  odd  with  odd. 

JEhere  are  other  selection  rules  which  are  not  v^id  in  general,  but  which  hold  to  the 
extent  that  the  intra-moleeular  interactions  can  be  represented  by  one  of  Miuid's  ideal 
coupling  cases.  \^en  each  state  participating  in  a  transition  belongs  either  to  ease  "a" 
Or  to  case  "b,"  then  A  ,  the  projected  electronic  orbital  angtdar  momentum,  obeys  the 
selection  rule  A  a  =  o  ,  +1  ,  -l  as  long  as  the  spin-orbit  and  rotation-electronic 
motion  interactions  are  small,  \^en  A  =  6  for  both  states,  there  is  the  further  rule 
that  the  sign  of  the  state  cannot  change*. 

2+**!;+  ^  r  ---2^  ,  but  not  2'^*— 2" 

Violations  of  the  above  rule  occur  when  the  rotation-electronic  motion  interaction 
cannot  be  regarded  as  small.  If  the  spin-orbit  interaction  is  smaU,  there  is  also  a 
selection  rtde  for  S  ,  namely  that  AS  0  . 


In  the  special  circumstance  that  both  the  tg>per  and  the  lower  state  belong  to  case  "a," 
the  following  two  selection  rules  hold: 

At  ^  0 

Afl  =  0  ,  +1,  -1  (but  not  n  ^  0  •=*  fi  =?  0) 

The  first  of  these  rules  holds  only  for  small  spin  orbit  interaction,  but  the  second  is 
not  so  restricted. 
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^  both  the  upper  and  the  lower  states  belong  to  case  "b,"  then  K  ,  the  quMtuna  number 
for  the  total  angular  momentum  apart  from  the  spin^  obeys  the  following  rule 

4K  =  b,  +1*  ^1  (but  not  K  =  0  K  =  0  for  S  S  ) 

4. 5  INTENSITY  FACTORS 


Tables  of  intensity  factors  (the  so-called  Mbnl- London  factors,  )  are  given  by 
Jevons  (Refi  8,  pp.  133“  139)  for  Hund's  coupling  cases  "a"  and  "b"  and  for  various  types 
of  transitions 4  For  intensity  factors  arising  from  other  coupling  conditions,  refer  to  the 
discussions  of  the  individual  molecular  systemSi  Table  A-2  lists  intensity  factors  for 
various  transitions  which  are  considered  in  the  calculation. 


According  to  the  rotational  sum  rule,  the  sums  of  the  line  strengths  of  all  the  trans¬ 
itions  from  or  to  a  given  rotational  level  are  proportional  to  the  statistical  weight  of 
that  level.  This  rule  affords  a  valuable  check  when  fine  structure  is  to  be  i^ored.  As 
an  example,  consider  a  11  «  T  transitions  Where  the  molecule  i$  C0Tig>led  according 
to  Hund's  case  "b."  Using  the  intcns%  factors  from  Table  lO(ii)  in  Jevons  (Ref.  8), 

We  combine  the  factors  for  the  lines  arising  from  level  J"  ^  K'  -  : 

RgfK  =  K"  +  1)  +  =  K")  +  QglK  ^  K")  +  -  K")  +  PgfK  =  K") 


K"tK"  +  2>  .  1  .  (K'\  +  1)  .  (K"  -  1)  .  (K"  -  1)^ 

"  2K"  +  1  2K"^  +  1  2K"  r  1  ^  ^  ^  -  1 


2K" 


which  is  equal  to  the  degeneracy  of  this  level,  i.e. ,  2  J"  +  1  =  2|k"  «  +  1  ^  21^'  . 

For  the  J"  =  K"  ^  level  we  obtain; 

Rj(K  -  K"  +  1)  ^  (K  +  K"  +  1)  ^  Q^(K  -  K")  +  (K  “  K")  +  Pj^(K  »  K**) 


_  (K"  +  2)^  .  (K"  +  2)  .  K"{2K"  +  3)  ,  1  _  .  (K"  +  IKK"  -  li 

■  2K"  +  3  (2K"  ♦  D^K''  +  3)  21?’  +  1  2K”  +  1  21?*  +  1 

-  2(K"  +  1) 
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which  is  equal  to  the  degeneracy  of  the  level,  viz.  ,  2J"  +  1 
Table  10  (ii)  in  Jevons  therefore  satisfies  the  sum  rule  and  is 
Other  results  may  be  checked  in  the  st^e  fashion. 

4. 6  THE  1  ^  -  X  (SCHUMANN-RUNGE)  SYSTEM  OF 

u  S  w 

4 .6.1  General  '  — — - 

With  Small  dispersion,  S-  X  bands  consist  Of  a  single  R  and  a  single  P  branch. 
Larger  dispersion  shows  each  line  to  be  resolved  into  three  components  of  about  the 
same  intensity.  There  are,  thus,  six  main  branches.  Six  very  weak  satellite  branches 
for  which  AJ  ^  aK  are  also  present.  Alternate  triplets  in  the  P  and  R  branches 
are  missing,  since  the  nuclear  Spin  of  the  oxygen  atom  is  zero.  Even  K  rotation 
levels  are  missing  in  the  X  state,  and  odd  K  levels  are  missing  in  the  B  state.  In 

our  calculation,  lines  will  be  treated  as  unresolved.  Figure  4-3  shows  two  groups  of 
-  R  R  P  P 

lines:  (  Rp  Rg,  Rg.  Qgi*  %2  ^  ^  ^1’  ^2’  ^3’  ^23 »  %2  ^  ‘  These  are  con¬ 

sidered  to  be  composite  lines  whose  Hdnl- London  factors  are  obtained  by  simply  adding 
\g)  the  Sj's  for  the  individual  Components  (Ref.  8,  pp.  180’- 181).  These  factors  turn 
out  to  be  3K"  for  the  P  branch  lines  and  3(K"  -i*  1)  for  the  R  branch  lines. 

4.6.2  Line  Frequencies 

3 

Rotational  terms  of  both  E  states  are  given  by 

F^(K)  ^  K(K  +  1)  -  K^CK  +  1)^ 

Table  A-9  lists  By  and  as  well  as  G^(v)  used  in  the  calculation.  The  included 
bands  and  their  corresponding  Franck-Condon  factors  are  given  in  Table  A- 10. 
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=  2|k’'  +  ij  +  1  =  2(K''  +  1)  . 
presumably  correct. 


4,6.3  Line  Mtensities 

For  tlie  Oo  Schumann-Runge  system,  then,  the  Hj„  .  jt's  6^^®’* 


16  7t^  Lq 

3  he  *'k'*.K‘ 


q(v’,V’') 


K"  +  1  if  K'  =  K”  +  1 
K««  if  K'  =  K"  -  1 


where  K*'  ~  i,  3,  5,  .  •  • 
and  Sj  -  3  . 

4,7  THE  B^n  (FIRST  POSITIVE)  SYSTEM  OF  Ng 

s 

The  structure  of  bands  is  very  complicated,  particularly  if,  as  is  the  case  here, 

the  coupling  in  the  t  state  is  near  Hund‘s  case  «a.«-  Since  there  are  three  sub^bands, 

4  ^  h*  %  -  .  %  -h*  ,  and  each  sub^band  has  nine  branches  (three  for 

each  triplet  component  of  the"lower  state),  there  are,  in  all,  27  branches.  Of  these 

we  expect  9  branches  to  be  relatively  strong  (  Pj,  ^2’  ^2’  ^2’  ^3’  *^3  ^  ’ 

ten  satellUe  branches  to  be  somewhat  weaker,  ai^  eight  to  be  weak  branches 
o  o^  o^  s„  Sp.  s_  np  Tr  ,  ).  We  will  omit  the  latter  ei|d>t 
(  Pl2'  %3'  ^23>  ^21'  %’  %’  -13’  ^31^- 

branches  entirely.  The  19  remaining  branches,  along  with  their  corresponding 
intensity  factors,  are  listed  in  Table  A''2.  (See  also  Fig,  4*4.) 

Formidas  for  tte  rotational  levels  of  a  ^2  state  have  been  derived  by  ScMapp  (Ref.  9): 
F  (K)  .  B  K(K*  1).  (2K.3)B  ♦  TW  ^  D 

Fv2(K)  ^  B^K(K+  1)  _ _ 

F^gfK)  *  K(K+  1)  -  (2K  *  1)  -  A  +  ^ {2K  ^  VT  -2X8^  -  Yk 
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where  ,  and  F^^  refer  to  the  levels  with  J  =  K  +  1  ,  K  ,  and  K  -  1 

respectivelyi  ahd  where  a  and  y  are  (small)  constants  ^  From  the  above  formulas 
must  be  snbtracted  the  usual  term  (K  +  1)^  to  allow  for  centrifugal  distortion. 

3  -  + 

Now  in  the  ease  Of  A  2  State  of  N„  ,  both  splitting  constants  are  very  small  andi 
furthermore,  are  known  only  very  approximately  for  most  of  the  vibrational  levels, 
"Therefore,  ^  areTuStmed  in  ignoring  the  spin  spiming  and  tisihg,  instead,  the  term 
formula  for  ^2  states: 

F^(K)  K(K  +  1)  -  K^(K  +  1)^ 

The  values  of  and  for  each  of  the  vibrational  levels  considered  are  given  in 
Table  A-3. 


Q 

Rotational  term  formulas  for  the  n  state  are  from  Budo  (Refs.  10,  11)  and  are  for 

'O 

any  d^ee  of  uncoupling: 


if 

=  ®v  [«•'  *  •  ®y('  *  if 

F,3(j)  •  [jW  +  1)  +  ^ 


where 

^  Y^(Yy  -  4)  +  3/4  +  4J(J  +  1) 
Zg  ^  Y^(Y^  -  1)  -  4/9  ^  2J(J  +  1) 


and  Y^  ^  ^  measure  of  the  degree  of  cotpling  of  the  spin  to  the  intemuclear 

axis.  For  large  rotation  ^vl-  F^2  '  ^y3  f*’  ^  series 
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with  J  =  K  +  1  ,  K  ,  and'  K  -  1  reapeGtiveiy.  Note  that  levels  ’ 

and  F^3(1)  <io  not  exist  beGause  of  the  requirement  that  J  s  n  =  (  a  +  2 1  for  Hund's 
ease  "a."  The  1  *  D  ,  and  Y  values  for  each  vibrational  level  are  given  in 
Table  A- 3.  Sinee  A  type  doubling  is  small  except  for  very  large  K  ,  it  has  been 
neglected  in  the  above  formulas^ 


The  vibrational  term  G  (v)  is  also  listed  in  Table  A-3 

o 

system  is  given  by 


Hj„  for  the  first  positive 


Lh 

he 


i/|  R  r  q(v',v")  x 
6 


2/3  if  K"  even 
4/3  if  K"  odd  , 


4. 8  THE  N2  C  ®  (SECOND  POSITIVE)  SYSTEM 
4.8.1  General 


The  structure  of  bands  is  simple  if  both  states  belong  either  to  Hund's  ease  "a" 

or  to ^Cmd's  case  "b."  If  both  states  belong  to  ease  "a,"  the  selection  rule  A2  -  0 
allows  division  into  three  sub-bands:  ,  and  •  Should  A  lype 

doubling  be  disregarded,  each  sub-band  has  a  strong  R  ,  a  strong  P  ,  and  (except  for 
^  ^  wealt  Q  branch.  If  both  %  states  belong  to  case  "b"  or  if  both  go  over 
from  case  "a"  to  case  "b"  with  increasing  rotation,  the  same  six  strong  bands  occur. 

For  aU  values  of  K  in  case  "a,"  and  for  large  values  of  K  hi  case  "b,”  the  three  P 
branches  are  close  together  and  the  three  R  branches  are  close  together,  giving  rise 
to  a  characteristic  triplet  structure. 


4,8,2  Line  Frequencies 

In  the  second  positive  system,  both  %  states  belong  to  case  "a"  for  small  K 
values.  Both  go  over  to  case  "b**  for  large  values  of  K  .  Figure  4-5  shows  the  six 
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3-3- 

main  braAckes  of  the  C  n  -1  n  system  and  the  Gorrespondtag  energy  level  diagram  i 

n  g 

A  rapid  transition  from  case  ''a”  to  Gase  "b"  has  been  assumed.  Note  that  both  states 

3-3-3 

are  "regular,"  i.e.  ^  n„  <'^n.  The  a -type  doubling  shown  in  the  figure  is 

O  Ji  ib 

greatly  exaggerated.  (  A^type  doubling  will  be  negleGted  in  this  study;  the  formulas 


for  the  rotational  terms  f , 


F.  ,  and  F«  should  not  be  affected  other  than 

2t  o 


trivially. )  A  few  of  the  levels  are  labeled  with  the  associated  nuclear  spin  statistical 
Weights  which  is  Written  to  the  right  of  the  level.  Since  nitrogen  nuclei  follow  Bose 
"Statistics  (t  -  1)  ,  the  symmetrical  levels  (s^^^ill  haVe^  the  Mgher  siatistical  weight 
ci>j(s)  =  (21  +  l)(l  +  1)  j  and  the  antisyMmetrical  levels  (a)  will  have  the  lower  statis¬ 
tical  weight  w*(a)  -  (2l  +  1)  I ;  i.  e. ,  o)  (s)  “  6  and  a;  Ja)  =  3  .  Therefore  ^e  total 

Weight  for  the  A  doublet  is  co^  =  9  ,  which  is  just  ^  Of  the  "maximum"  weight 

2  ~  * 

0,^21  H  ^  1)  =  is  that  one  wotdd  e:!^et  for  a  heteronuclear  system  of  the  same  ^ipe 

and  for  which  both  nuclei  have  ^  ^  ^  * 

The  lines  indicated  in  the  figure  are  in  agreement  with  the  selection  rules 
S  s  ,  a  —  a  ,  and  J  -  o  ,  *  i  ;  except  that  Ad  -  0  is  forbidden  for  n  -  0  —  fl  =  6  . 
(There  can  be  no  Q  branch  for  n^-  fl-  ,)  The  Qg  and  branches  have  line 
intensities  that  fall  off  as  1/ J"  ,  and  hence  the  total  intensity  residing  in  these 
branches  is  Small.  Therefore,  these  branches  will  be  omitted  in  the  present  study. 
Finally,  we  note  that  the  line  F,(o)  is  missing  in  n  -  fl  ,  the  lines  P„(o)  , 

^^(l)  ,  and  R^{o)  in  ,  and  the  Itoes  Pglo)  ,  Pg(l)  ,  Pg(2)  ,  Rg(o)  ,  and 

R  (1)  in  in  conformity  with  the  rule  that  J  s  n . 

V  « 

Formulas  for  the  rotational  terms  ,  Fg  ,  and  Fg  have  been  given  by  Budo^ 

(Ref.  10); 


Fj(J) 


FgCJ) 


FgfJ) 


J(J  ♦  1)  +  ^  -  I  ^1  -  B^(J  ♦  3/2)*  (J 
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where 


^  Y^(Y^-4)  +  |+4J(J  +  1) 

z.  ^  AY JY^  -  1)  “  I  “  2J(J  +  1)  , 


and  Y  »  A^/  is  a  measure  of  the  strength  of  cot^ling  of  the  spin  to  the  intemuclear 

axiSi 


If  ^y«i  yi  =  ^oo  +  “  G|J(v*’)  ,  then  the  frequencies  of  the  lines  will  be  given  by 


V  "  V.v.j-(J’-:Hi)  "  V.v  *  *  D 


p 

V  ^  =  V 
*^.T"  v"  V* 


J"  “  1)  V’.V  h 


^  *  n(J"  "  i)  -  W) 


for  n  ^  l«  2,  3  .  Again,  A  type  doubling  is  neglected. 

4.8.3  LUie  Intensities 

The  two  lines  arising  from  the  components  of  a  A  type  doublet  overlap  in  frequency 
and  are  identical  except  for  the  ,  one  having  the  value  6  (symmetric)  and  the  other 
3  (antisymmetric).  This  allows  us  to  construct  the  H „  ,  „  ,  ~  tor 

the  whole  doublet  by  simply  adding  the  H's  for  the  individual  components,  viz; 


Sjt-  L 
3hc  2 


^  *  v.j’  i^e'^  ®j"  [m  *  m] 


1-0  2 

“  She  >',3'  *  <lty',V’)  Sj,, 


where  9/2  is  the  value  of  Q 
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.  3—  3 - 

UeM- London  factors  for  n-  n  systewis  under  various  coupling  conditions  have  been 
given  by  ludO'(Ref.  12).  For  this  study  it  will  be  sufficiently  accurate  to  use  the 
Hbid.- London  factors  for  a  tj^e  transition: 


5^ 

’ji. 


=  jw  1 


=  J**  . 


The  final  formulas  for  v’.jV.  J”,  J*  then  become 


R.  8ff  -  L^  R^  » 

"j.?  iRg* 


and 


P-  L  P  o 

^  ‘'j.“  I  Re* 


for  n  e  1,  2,  3  .  FlnaRy, 


V.j"  =  *k  [®o<’">  ♦  ■>  -  1.  2.  3 


4. 8. 4  MolectRar  Constants 


3  3 

The  molecular  constants  for  the  C  n  and  B  n  states  of  N„  which  were  used  in 

u  g  2  ^ 

this  study  are  given  in  Table  A-3.  The  values  of  ,  and  are  due  to  Budo. 

The  vibrational  terms  G^(v)  were  taken  from  Coster,  Brans,  and  Van  der  Zlel 


(Ref,  13), 


The  quantity  F_  listed  to  Table  A-3  gives  the  maximum  value  that  any  of  the 
rotational  terms  F  (n  ^  1,  2,  3)  may  have.  The  rotational  formulas  are  probably 

Vjll 

not  accurate  beyond  J  =  50  or  60  .  The  extent  to  uddch  they  fail  is  not  known,  but 
it  is  felt  not  to  be  great  until  about  J  =  80  . 
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fable  A‘^4b  shows  the  bands  included  iti  the  present  calculation  along  with  the  Franck^ 
Condon  factors  usedi  Note  that  the  vibrational  levels  in  the  upper  state  stop  at 
V'  =  4  i  The  potential  curve  does  not  break  off  at  v'  =  4  however,  and  although 
transitions  from  v'  >  4  have  not  been  observed,  one  shotdd  in  a  local  thermodynamic 
equilibrium  situation  go  to  v*  >  4  ^  Franck^^Condon  factors  and  S'  values  are  not 
Imown  for  v'  >  4  ,  but  me  S's  could  be  obtained  by  extrapolation, 

-4:9  THE  B  h*  -  X  (Fmaf  negative)  system  of  ^ 

4.9.1  General 
2 

Since  £  stMes  always  belong  strictly  to  Hund's  coupliiig  case  giving  rise  to  the 

-  -  -  22- 

selection  ride  AK  -  ±  1  with  AK  s  0  being  forbidden,  the  structure  of  a  S 

transition  is  Simple.  There  are  two  prMclpal  branches,  an  R  branch  mid  a  P  branch 
each  of  which  may  be  resolved  into  three  components  according  to  the  rule 
AJ  =  0  ,  ±  1  ,  For  one  of  these  (AJ  =  0)  AJ  AK  . 

4. 9. 2  Line  Frequencies 

The  energy  level  diagram  and  transitiims  which  iq>ply  to  first  negative  bands  are  snoi^ 
jtai  Fig.  4-6.  Rotsttan  terms  of  lower  and  tq;>per  E^  states  are  both  given 

F^(K)  =  S^  K(K  +  1)  -  k\k  +  1)2 

where  the  spin  splitthig  of  each  rotational  level  is  ignored.  Table  A-5  lists  the  values 
of  ,  and  G^(v)  used  in  the  calculation.  The  included  bands  and  corres¬ 

ponding  Franck-Condon  factors  ture  given  in  Table  A-6.  Bands  with  very  small  F-C 
i^ictors  have  been  omitted. 

The  Unes  ud  tiieir  correspoiding  Hbnl- London  factors  are  given  in  Table  A-2. 
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Fig.  4->6  B  Sj*  X  Band  of  Ng  (1st  Negative  ^sten^ 
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4.9^3  Ltoe  Mtensities 


2-57-62“2 


Grouping  together  ,  and  into  a  single  line,  the  Sum  of  the  Sj's  is 

2K''  .  For  the  group  ,  R^  ,  and  S’*®  Strength  faetors  is 

2(K"  ^  1)  ^  goring  Ihe  spin,  the  rotational  levels  can  be  labeled  by  K  .  The 


H 


's  are: 


167f^  L 
<3 

She 


V.K  '*e'  *»•  ■ 


1/3  if  Kl'  odd 


adtere 

=  K"  +  1  if  K*  =  K**  +  1 
=  K?*  if  K'  =  K**  “  1 

The  factors  to  allow  for  K  odd  or  even  arise  because  of  the  effect  of  syMmetry 
properties  on  t-t  transitions  of  hoMonuclear  molecules.  Even  K"  levels  are 
symmetric  (Wj(s)  »  6) ,  vdiile  odd  K"  levels  are  antisymmetric  (w^fs)  «  3). 

4.10  THE  B^n-X^n  (BETA)  SYSTEM  OF  NO 

4. 10, 1  Genera 

2  2 

The  structure  of  a  n«  n  system  is  simple  if  both  states  belong  to  Hand's  coupling 

case  "a"  or  both  belong  to  case  ”b''  or  if,  as  to  the  present  case,  both  states  are 

totermediate  between  cases  "a”  and  "b”  and  go  over  together  from  case  "a"  to  case  "b” 

with  increasing  nuclear  rotation.  There  are  two  sub^bands  to  the  NO  Beta  system: 

2  2  2  2 

^1/2'”  ®1/ 2  ^3/2’ ®3/ 2'  Eu®h  sub-band  consists  of  a  strong  R,  a  strong 

p  ,  and  weak  Q  branch.*  Both  the  B  and  X  states  are  regular,  that  is, 

2  2  —1 

^1/2  ^3/  2  '  intensity  of  the  Q  bruidies  falto  off  like  J  ,  we  will 

neglect  these  branches  entirely. 


♦See  Fig.  4-7. 
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Fig.  4-7  B 


-  X  Band  of  NO  {p  System) 
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4. 10. 2  Like  Frequencies 


2-  -  2 

Both  the  B  II  and  X  n  states  of  NO  are  intermediate  between  Hund's  coupling 

cases  "a"  and  "bj"  with  the  X  state  being  nearer  case  "a"  (A  -  124. 2  cm  )  than 

the  B  state  (for  which  A  ~  32  cm  ^  ).  A  general  formt^a  for  the  term  values  of  a 
2— 

n  state  with  coupling  intermediate  between  cases  "a"  and  "b**  has  been  derived  by 
HUl  and  Van  ^^eck  (Ref.  14): 


T^Cv.J) 


te  +  G^(v)  +  B^  (j  ^  1  +  (-  I)*"  KJ)  - 


(n  ^  1) 


(J  +  1)^  (n  ^  2) 


where 


2 

-’f,  '^V  * 


2  -  2 

Here  the  subscript  n  =  1  stands  for  a  and  n  =  2  stands  for  113^2  ‘ 

subscript  n  has  been  affixed  to  the  vibrational  term  G(v)  ,  since  the  different  suh^" 
states  of  the  multiplet  may  have  sii^Uy  different  vibrational  levels.  Formulas  for 
the  line  frequeneies  of  the  two  R  and  two  P  branches  follow. 


Hl(J)  ^  »'e  ♦  +  F’^(J"  +  1)  -  F'^(J") 

R2(J)  ^  »'e  +  »'v  ^  +  1)  •* 


Pl(^ 


=  |/e  v\  ’  +  F'j(r  -  1)  -  P{(J") 


uiiere 


VJ.J)  =  +  F^(J"  -  1)  -  F”(J”) 


Tt  «  'pt 

e  e 
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and 

yj)  ^  F^(V,J)  1^ 

» 

(J  +  TT  n“^"2 

k 

with  1^  ,  D^  i  G^(v)  ,  and  G2(v)  being  given  in  Table  A-7. 
4il0.3  Line  Intensities 


t  ^ 

1/2 

V  *  5/  •  1  ♦  »>“ 

hi) 

2- 

As  the  n  states  approach  case  "b,"  the  Selection  rule  AK  =  0*41  holds.  Also 
branches  with  AK  ^  AJ  are  very  wes^.  Disregarding  A  type  doubling  and  neglecting 
the  weak  Q  branches  and  satellite  branches,  the  band  structure  is  siniilar  to  that  of 
a  2-  2  transition,  i.  e. ,  4  strong  branches,  intensity  factors  for  case  "a"  and  case 
••b"  couplhig  are  given  in  Table  A-2.  If  the  satellite  lines  are  considered  as  part  of  a 
composite  lines,  but  Q  branches  are  still  ignored,  the  intensity  factors  may  be  given 
with  sufficient  accuracy  for  our  purpose  by: 


R, 


-J" 


O'*  +  1 


J” 


and  the  H  factors  are  then 


The  Franck-Condon  factors  q(v',v'')  are  presented  in  Table  A-8a. 
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THE  A  -  X  (GAMMA)  SYSTEM  OF  NO 

4. 11. 1  General 

_  2  _  _ 

The  X  n  state  of  NO  belongs  neither  strictly  to  coupling  case  "a"  nor  to  case  "b," 

but  to  a  transition  case  which  goes  from  case  "a"  to  case  "b"  wiih  increasing  rotation. 

2 

However j  the  A  s  state  belongs  strictly  to  case  "b,"  and  as  a  result  of  this  com- 

.  _  -  2 

tiinatlon  Me  structure  of  ^e  NO  gamma  bands  is  quite  complioatod. — The — -H-  state  has 

been  discussed  previously  (see  NO  Beta  system).  Figure  4-8  shows  levels  and 
transitions  applicable  to  Me  NO  gamma  System. 

4.11.2  LMe  Frequencies 

The  general  term  formula  is  given  in  the  preceding  discussion  of  the  NO  Beta  system. 

j  ,  G^^\v)  ,  and  Gq^^v)  for  the  X  state  are  presented  in  Table  A-7  ^ 

as  are  B  ,  D.  ,  and  G  (V)  for  the  A  State.  The  quantum  number  K  has  been 

V  V  O  Q 

formally  extended  to  levels  having  Small  rotation  in  the  ease  ^  the  A  £  state  (Bel.  5). 
Rotatiomal  terms  for  this  state  are  given  by: 

F^(K)  «  B^K(K  +  1)  -  D^k\k  +  1)2 


4. 11. 3  Line  Intensities 

Mtensity  factors  for  the  twelve  br^ches  are  given  (apart  from  an  arbitrary  constant) 
by  Earls  (Ref.  15).  These  formidas  were  normalized  to  obey  the  sum  rule.  The 
resulting  expressions  are  given  in  Table  A”2.  After  cpmbMMg  brMiches  and  satellites 
and  omitting  two  weak  transitions,  6  composite  lines  serve  the  needs  of  our  calculatioii. 
btensity  factors  for  these  lines  are  also  given  in  Table  A-2. 
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where  the  are  the  composite  toteasity  factors  described  above  <  The  Franck- 
CoiMiop  factors  q(v'  .v*')  for  the  transitions  are  listed  in  Table  A-8b. 
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Section  5 

THE  SACHA*  CODES 

5. 1  MAGNETIC^TAPE  LINE  ATLASES 

Digital  GOffiputer  codes  (for  the  IBM  7090)  have  been  constructed  to  write  the 
Kflagnetie^tape  line  atlases  for  the  six  Molecular  systenrs  to  be  included  in  the  ealeu= 
lation>  For  each  system,  the  Sipectral  lines  and  their  "H"  and  "E"  functions  £see 
Eqs*  (3i2l)  and  (3. 22)  and  Subsections  4.6  to  4.  llj  are  generated  according  to  vibra^ 
tional  and  rotational  quantuM  nuMbers  and  labeled  by  an  identification  vector »  These 
lines  are  then  merged  and  sorted  accorduag  to  line  frequency,  For  each  spectrai 
line,  there  are  three  decimal  numbers  and  one  octal  number  (the  identification  vector) 
on  tape: 

p  a  E 

a  a 

(frequency)  (identification  vector) 

where  o  is  formed  of  the  upper  and  lower  vibrational  and  rotational  quantum  num-« 
bers,  the  branch  number,  and  an  integer  from  1  to  6  which  specifies  the  molecular 
system  to  which  the  line  belongs. 

The  systems,  their  "key"  integer  n ,  and  the  total  number  of  lines  stored  on  tape  for 
each  system  are; 


1?  -  1 

^2 

Schumann-Runge 

13,836 

2 

First  Positive 

58,476 

3 

^^2 

Second  Positive 

16,750 

4 

First  Negative 

21,306 

5 

NO 

Beta 

15,760 

6 

NO 

Gamma 

25.400 

Total  Lines 

151,528 

*l^ectral  Absorption  Coefficient  of  Heated  Air. 
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In  addition  to  the  individual  system  line  atlases,  a  tape  consisting  of  all  the  Systems 
merged  according  to  frequency  has  been  written.  This  merged  tape  containing  all 
151,528  lines  will  serve  as  the  spectral  line  atlas  for  air.  By  masking  techniques,  the 
origin  of  each  line  may  be  found  from  its  identification  vector,  i.e. ,  its  species,  lower 
level,  etc.  To  compute  the  absorption  coefficient,  the  atlas  tapes  must  be  used  in  eon- 
junction  with  population  numbers  and  partition  fiiinctions  correspondiag  to  the  desired 
air  temperature  and  deasity.  Population  factors  may  be  obtained  by  interpolating  in 
^4ablei^^^Gilmnr«4Be£-JLfi)  Thp  vihratinin^rntatinri  partition  funCtiona  are  giveil^ 
with  sufficient  accuracy  for  our  problem,  in  Table  A-lc^ 

All  of  the  Franck-Coadon  factors  used  in  the  atlas  construction  are  from  R.  W. 

Nlcholls  (Refs,  17^19).  The  electronic  f^^values  Were  taken  from  Treanor  and  Wurster 
(Ref,  20^Oq  Schumann- Runge)i  Bennett  and  Dalby  (Ref.  21 First  Negative)  i  and 

M  A 

Keck  et  al,  (Ref,  21-reniaining  band  systems). 

5.2  AN  AVERAGE  TRANSMISSION  CALCULATION 

As  a  first  effort  to  assay  the  value  of  the  method,  a  code  was  written  to  calculate  the 
average  transmission  of  optical  radiation  throng  an  isothermal,  homogeneous  slab  of 
pure  molecular  oxygen  as  a  function  of  distance  through  the  slab.  Variable  parameters 
in  the  calculation  are  the  slab  temperature  T  ,  the  line  half-width  o  ,  basic  fre¬ 
quency  intervals,  and  the  number  of  frequency  points  on  the  line  wings.  The  average 
transmission  at  average  frequency  v  is  given  by 

Av  L  0! 

where  ^he  contribution  per  particle  to  the  spectral  absorption  coefficient  at 

frequency  v  due  to  spectral  line  a ,  N  is  the  particle  density  of  molecules  in  the 

lower  electronic  state  for  the  Scbumann-Runge  transistions,  and  x  is  the  distance 

2  “3 

through  the  slab.  The  units  of  are  cm  ,  those  of  N  are  cm  ,  and  those 

of  X  are  cm.  T  is  the  center  frequency  of  the  interval  Av ,  where  v  has  units 
cm  When  expressed  in  these  imits,  the  frequency  is  generally  referred  to  as  the 
"wave  number. "  The  sum  runs  over  all  lines  a  u^ich  contribute  appreciably  at  v. 


j  dv  (5. 1) 
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The  cakulation  was  carried  out  by  computing  the  average  transmission  as  the  product 
Nx  was  variedi  By  using  this  technique ,  the  need  to  explicitly  Specify  N  Was  circurn'^ 
vented  for  these  preliminary  studies.  The  product  Nx  is  referred  to  herein  as  the 
"effective  depth. "  Since  these  preliminary  calculations  were  done  for  ten^eratures 
at  or  below  5,606'’K,  emission  was  neglected. 

Some  representative  results  of  the  calculation  are  presented  in  figs.  5"“!  through 
Figures  5°^!  and  5^2  demonstrate  the  effect  of  variation  of  the  line^nddth  parameter 
for  fixed  values  of  the  other  parameters  as  noted  in  the  figures. 

Average  transmission  is  plotted  versus  frequency  in  figs.  5^3  through  5^5,  for  various 
Values  of  T,a,  and  Nx.  Although  smooth  curves  have  been  drawn  ^  the  calculation 
was  carried  out  only  for  isolated  frequency  interv^s  of  100  cm  ^  widths  These 
intervals  were  centered  at  ten  frequency  points  from  31,500  em~^  to  49,  500  cm""^,  the 
difference  between  points  being  2,000  cm^^. 


Mien  only  (me  molecular  species  is  involved,  die  relatitm  between  particle  and  linear 
absorption  coefficients  is 


M(V)  =  Nic(v) 


(5.2) 


x^ere  p(v)  has  the  usual  units  cm 
5.3  COMPARISON  OF  RESULTS 

An  earlier  calculation  carried  out  at  this  laboratory  by  Meyerott,  Sokoloff,  and  Nicholls 
(Ref,  1)  affords  the  only  known  results  wMch  are  comparable  with  those  of  the  present 
work,  hi  the  earlier  calculation  (hereafter  referred  to  as  the  MSN  paper) ,  the  average 
(Uscrete  absorption  coefficient  was  computed  over  a  much  wider  frequency  range,  using 
basic  averaging  intervals  of  0.25  ev  (roughly  2,000  cm^^).  Computation  was  carried 
out  by  desk  calcidator  and  hicluded  a  large  number  of  vibrational  bands  for  each  of  the 
six  molecular  transitions  needed  for  the  air  absorption  calculation,  Rotational 
stiucture  was  not  included,  with  the  total  band  contribution  taken  at  band  head  energy. 
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Jig.  5-2  Effect  of  HaU-Width  Variation  cm  Optical  Properties  of  for  v  =  47500  Cm  and  T  =  2000“K 


I 

2-57-62-2 

I 

I  Since  this  work  was  published  previous  to  the  availability  of  many  of  the  Og 

Schumann-Runge  Franck-Condon  factors,  these  bad  to  be  estimated. 


Figure  5-6  was  prepared  by  ealouiating  average  transmissions,  utilizing  the  O 

2 

Schumann-Runge  absorption  coefficients  listed  in  the  MSN  paper,  after  modification 
of  the  electronic  f-number  and  level  population  to  agree  with  those  assumed  in  the 
present  work. 


The  cause  for  the  differences  between  the  two  Galculations  has  not  been  extensively 
investigated.  Apart  from  the  large  difference  between  basic  frequency  intervals 
employed  in  the  averaging  process,  the  major  contribution  to  the  discrepancy  appears 
to  arise  from  the  omission  of  several  important  band  systems  in  the  previous  work. 

Inclusion  of  many  more  bands,  consideration  of  rotational  structure,  and  smaller 
frequency  intervals  for  averaging  are  among  the  factors  giving  rise  to  the  more 
smoothly  varying  transmission  obtained  in  the  ^CHA  results.  More  Frsmck* 

Condon  factors  were  available  for  the  present  study,  and  better  values  had  replaced 
some  of  those  used  in  the  earlier  caiculations. 

5.4  A  CALCULATION  FOR  MEATH?  AIR 


I 

I 

I 


Another  machine  program  in  the  SACHA  series  is  nearing  completion.  This  code 
will  compute  the  average  transmission  of  optical  radiation  through  a  Slab  of  heated 
air.  Mean  absorption  coefficients  will  also  be  obtained  over  the  frequency  rmige  of 
interest.  The  six  molecular  transitions  discussed  in  this  report  vdll  be  included  in 
the  calculation,  with  the  species  concentration  and  fractional  population  for  the  lower 
electronic  state  of  the  transition  being  obtained  by  interpolation  in  Gilmore's  Tables. 
At  a  later  date  provision  may  also  be  made  to  include  the  continuum  contribution  to 
the  absorption  coefficient.  Transmissions.and  absorption  coefficients  will  also  be 
obtained  for  the  individual  constituents  of  air. 


I 

I 

I 
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f^le  A-1 

MOLECULAR  SYSTEMS  CONSTANTS  (GENERAL) 

a 

Spectroscopic  Qu^tities 


System 

Transition 

— rxcm  ~poi — 

f 

“  IR  T"  . . 

(Absorption) 

f  Measurement 

I«gl 

O  Schl>*nann— RUnge 

X®2"“-B®2"‘ 

29,300 

O.048 

■  36 

3,  293  X  10 

2 

No  First  Positive 

2 

g  u 

aX-b% 

12, 500 

0,02 

3.34  X  10"^® 

Nq  Second  Positive 

2 

b\-cX 

32,000 

0.07 

4. 6  X  lo"®® 

N^  First  NegaUve 

2 

25,500 

0,0348 

2.85  X  lO"®® 

NO  Beta 

X  ^  1 

23,800 

0,008 

0,7  X  lO"®® 

NO  Gamma 

xS  — A  H 

38,500 

0,0025 

0. 14  X  10‘®® 

b 

Statistical  Quantities 

Molecule  and  State  ^ 

^S 

‘"a 

o 

0 

0  2 

3 

1 

2  g 

N  A  '2* 

1 

1  2 

3 

1 

2  U 

N_  B  'H 

1 

1  2 

3 

2 

2  g 

N*  X  -2* 

1 

1  2 

2 

1 

-  2  g 

NO  x^n 

1 

0  1 

2 

2 

A-2 

LOCKHEED  MISSILES  ft  SPACE  COMPANY 
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Partition  Functions 

MoleGule  and  State  Vibration-Rotation  Partition  Functions 


02  **2:; 

Ql(T) 

0.4834  T(1  +  0. 0000149  f) 

1  -  exp  {-  2256/f ) 

%(T) 

0.48468  Til  +  0.  0000174  T). 
“  1  -  exp  (-  2081. 2/T) 

**2 » \ 

QgCT) 

0.42670  Til  +  0.0000144  t) 
1  -  exp  (-  2474/T) 

Q4(T) 

0.3617  T(1  +  0.0000102  T) 

1  -  exp  (-  3152.  5/T) 

NO  X^H 

Qg(T) 

0.4698  Til  +  0.0000119  T) 

1  -  exp  (“  2719/T) 
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Qi(J) 


Q 


Ri2(J‘1) 


Pi(J  +  1) 


Qi2<^) 


Ri3(J^1) 


1) 


[(r  +  j  -  i)u^  +  (Y  -  2)  *  2j(r  ^  i)j 


JC^(J) 


<j^-l)[(J^^  l)(Y^2)-uJ 


JC^(J) 


j[j(J  +  2)Uj  +  (J  +  2)  (Y  “  2)  +2  (J  -  1)  (J  +  1) 

(j  + 1)  ■r3>c^<jy 


[(J 


(2J  +  1)  (J^+  J  ^  1)  (Y  ^  2)  + 


j(j  +  i)e^(j) 

(J  ^  1)^  (J  +  1)[(J  +  l)u^  "  (Y  ‘  2)  2J(J  +  1) 


J(2J-i)C^(J) 


RoiJ'i) 

8(J  -  1)  (J  +  ir 

A 

JCgTJ) 

2[j(Y  -  2)  -  2]^ 

JCgR) 

Q2(J) 

8(2J  +  1)  (J^  +  J  - 

2jC<J^^  H  Y^  2  123+  3)3" 

(j  +  i)(2j  +  ayCgitJ) 
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Transition  and  Branoh 


Table  A*2  (oontd) 


2(J  -f  lvLj(Y  *41+2. 
J{2J*  1)G2(J) 


Table  A«2  (coat'd) 


rransitlQn  and  Branch 


where 

u. 


u* 


[y(Y  -  4)  +  4J^]  ^ 
[y(Y-4)  +  4(J+ 


C^(j)  =  J(J+  1)Y(Y  -4)  +  2(2J+  1)(J  -  1)J(J+  1) 
e2(J)  ^  Y(Y  -4)  +  4J(J  +  1) 

C3(J)  =  CJ  *  l)tJ  +  2)Y(Y  -  4)  +2(2J  +  1)J(J  +  1)(J  +  2) 


P^(J  =  K+  1/2) 

( J  =  K  +  1/2) 
®Q2j(J  *K  +  1/2) 


P2(J=K^1/2) 


-  1/4  ^  2K(K+  IV 
J  2K+  1 

1)^  -  1/4  ^  2(K+  1)  (Kh 
J V 1  ■  2  K  +  $ 

2J  *  1  _ 2  (K+  11 

4J(J  +  i)  ”  (2K+  l)(2K+  3) 

^  1/4  ^  2K(  K  -  1) 

J  ^  2  K  -  1 


R2(J=^K-1/2) 

*’q^2(J  ^  K  ^  1/2) 
^ere 


<J^  11^^  1/4  2K(K^  1> 

J+i  ~  2K+1 

2 J  4  1  2K 

4j  (J  +  1)  “  4  RZ  -  1 


Pj(K  +  1/2)  +  PgCK  -  1/2)  +  (K  -  1/2)  -  2K 


R,(K  +  1/2)  +  R-(K  -  1/2)  +  Rq  (K  +  1/2)  =  2(K  +  1) 
1  Z  '^^l 

A=6 
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Table  A-2  (cont'd) 


Transition  and  Branch 

2n  *  2jj 


Ri„(J)  orR_(j) 


laJJ) 


Qop(J)  or%,(J) 


P2e(J)  or  P„  .<  J) 


%,(J) 


Case  "a"  -  Case  "b” 


Case  ''b"  -  Case  "b” 


J(J  +  2) 
J  ♦  1 


i2J  +  n 
J(J  +  1) 


(J^  -  1) 
j 


2K(K+^ 


(FTTH^K  +  3)  K  -  J  -  1/2 


2K^(K  4-  2) 

(K  +  1)(2K  +  T) 

_  2(2K  +  31„ 


(K+  1)(2K  +  1) 


2(2K  -  1) 
K(2K  +  1) 


-  K  -j-  ir/Z 


where  K  =  j  -  1/2 


where  K  =  J  +  1/2 


2(K  -t-  ^  n  where  K  =  J  -  1/2 

K(2k:  +  i) 


2(K  ^  1)^(K  »  1) 
K(2k  -  1) 

2(K^  1) 


where  K  =  J  +  1/2 


where  K  =  J  1/2 


2K<K  +  2) 

(k  +  l)  (2ic +  T)  (2k  +1)  where  K  *  j  -  1/2 
2 

K(k  +  i)(2k  +  1)  ^^®re  K  s  J  -  1/2 


1)(2K  +  1)  K  ^  J  +  1/2 
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Table  A-2  (coat'd) 


(2J  ♦  1)^  *  (2J  -I-  1)  U  (4J^  »  4J  1  »  2Y) 
8(J  +  1) 


Rl(J) 

QoCJ) 


Qi(J) 

P2(J) 

%<•» 

Pl2f» 

Pl(J) 


where 


Then 


:2J  1)^  (2J  1)  U  (45^ 

8(J  +  1) 


(2J  »  1) 


♦  4j  -  1)  ±  U(8J^  +  12J^  ^  2J  +  1  -  2Y 


8J(J  +  1) 


(2J  »  1)1(4J^  »  4 j  -  1)  »  11(8/  »  12J^  2J  *  7  ■!•  2Y 

8J(J  +  1) 

(2J  ♦  1)2  *  (2J  ♦  1)  U  (4J^  +  4J  -  7  2Y) 

8J 


(2J  »  1)^  ^  (2J  »  1)  U  (4J^  »  4J  »  1  -  2Y) 

8J 


U  =  -  4Y  +  (2J  + 


RjlJ)  +  %i(J)  -  ^  {ej  ^  I  ^  y(4J^  +  4J  +  1  -  2Y)} 
Qj(J)  +  ^P2l(J)  ^  jej  +  7  +  U(4J^  +  4J  +  1  -  2Y)J 


2J  +  1 

1 

8J 

1 

2J  + 

1 

8(J  + 

1) 

2J  +  1 

{ 

^  8J 

\ 

2J  * 

1 

8(  J  ♦ 

1) 

6J  -  1  ♦  U(4r  +  4J  +  1  *  2Y) 


Pj(J)  »  ^(j)  >  S  vn  (w  »  7  -  0(4J^  ♦  «  »  1  -  2T)} 
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Table  A-3 


State 


N2  MOLECULAR  CONSTANTS 


V 

By(Gm"S 

D^(10“®  Gin^S 

F 

max 

Go(v)(em'^) 

0 

1.8154 

6.0 

21.5 

16579 

0.0 

1 

1.7932 

6.0 

21.5 

14585 

1994. 2 

2 

1.7682 

6.0 

21.4 

12645 

3934. 4 

3 

1.7407 

7.5 

21.1 

10771.5 

5808. 1 

4 

1.7012 

11.0 

20.3 

8989. 8 

7589.8 

0 

1.6285 

5.8 

25. 9 

2  X  10“^ 

0.0 

1 

1.6108 

6.0 

26.2 

2  X  10“* 

1705.2 

2 

1.5925 

6.0 

26.4 

2  X  10“^ 

3381.4 

3 

1.  5735 

6.0 

26.8 

2  X  lO'* 

5028. 0 

4 

1.  5554 

6.0 

27.0 

2  X  10“^ 

6646. 6 

S 

1. 5364 

6.0 

27.3 

2  X  10'^ 

8236. 0 

6 

1.5172 

6.0 

27.6 

2  X  lO"^ 

9796,9 

7 

1. 4954 

6.0 

27.95 

2  X  10^ 

11328.4 

8 

1.4765 

6.0 

28.25 

2  X  lo"^ 

12831.0 

9 

1.4576 

6.0 

28.5 

2  X  10^ 

14304.7 

10 

1.4387 

6.0 

28.8 

2  X  10^ 

15749. 4 
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fable  A-3 

(contlniied) 

D^(l0-®cm'S 

Gq(v)  (em"^) 

1. 433 

5.5 

0.0 

100 

1.421 

5.5 

1432. 5 

100 

1.408 

5.5 

2837.0 

100 

1.395 

5.5 

4213.3 

100 

1.382 

5.5 

5561.4 

100 

1. 369 

5.5 

6880.9 

loo 

1.356 

5.5 

8171.9 

lOO 

1.343 

5.5 

9434. 1 

100 

1.330 

5.5 

10667.4 

100 

1.317 

5.5 

11871.7 

100 

1.304 

5.5 

13046. 7 

100 

1.291 

5.5 

14192.4 

100 
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FRANCK-GONDON  FACTORS  FOR  THE  A  *  -  B 

u  g 

(FIRST  POSITIVE)  SYSTEM  OF  N* 

2 


II 

S. 

v 

V" 

JL 

6 

0 

.3382 

4 

1 

.1318 

0 

1 

.3248 

4 

2 

.  2738 

0 

2 

.1899 

4 

4 

.1139 

0 

3 

,  8857''“ 

4 

5 

.  8823 

1 

0 

.4064 

4 

8 

. 6361"^ 

1 

2 

.1032 

4 

9 

. 8349*^ 

1 

3 

.  1782 

4 

10 

.7511*^ 

1 

4 

.1450 

4 

11 

. 5510*^ 

1 

5 

. 8647^^ 

5 

2 

.  2106 

2 

0 

.1975 

5 

3 

.1808 

2 

1 

.2120 

5 

6 

.  1057 

2 

2 

.1132 

6 

2 

.6148*^ 

2 

4 

.7724^^ 

6 

3 

.  2605 

2 

5 

.1275 

6 

4 

. 8305" - 

2 

6 

.1127 

6 

5 

.1040 

2 

7 

.7496"^ 

6 

7 

.8078"- 

3 

0 

. 5014'^ 

6 

8 

. 6967*- 

3 

1 

.2987 

7 

3 

.1065 

3 

3 

.1623 

7 

4 

.2706 

3 

6 

. 6910'^ 

7 

6 

.1290 

3 

7 

.1009 

7 

9 

.7733*- 

3 

8 

.9079’^ 

8 

4 

.  1561 

3 

9 

.6406*- 

8 

5 

.2438 

8 

7 

.1158 

*  indicate  the  power  of  ten  to  which  the  fraction  is  to  be  reised,  throughout 
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table  A-4a  (Cont'd) 


y' 

v" 

9 

5 

.  2029 

9 

6 

.1918 

9 

8 

. 7862^^ 

9 

9 

.6773‘i 

16 

5 

.5569“i 

10 

6 

.2402 

10 

7 

.1300 

10 

8 

.5219'"^ 

10 

10 

.  8370“^ 

11 

6 

. 8265^^ 

11 

7 

.i630 

M 

8 

. 7?06"^ 

11 

9 

. 8779^1 

11 

11 

. 7897-1 

12 

7 

.1140 

12 

8 

.2688 

12 

10 

.  1109 

13 

8 

.1480 

13 

9 

.2575 
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Table  A-4b 

FRANCK-CONDON  FACTORS  FOR  -  B  ^11^ 
(SECOND  POSITIVE)  SYSTEM  OF  Ng 


v' 

V" 

q 

v' 

v” 

q 

0 

0 

.4493 

2 

7 

.3860"^ 

0 

1 

.3287 

3 

0 

.2363''^ 

0 

2 

.  1469 

3 

1 

.2515 

0 

3 

.5226"^ 

3 

2 

.1630 

1 

0 

.3899 

3 

3 

.1181 

1 

1 

.1868“^ 

3 

5 

.8891"^ 

1 

2 

.2038 

3 

6 

.1345 

1 

3 

.2003 

3 

7 

.  1062 

1 

4 

.1124 

3 

8 

.616r^ 

1 

5 

.4839^^ 

4 

1 

.6957^" 

2 

6 

.1349 

4 

2 

.3034 

2 

1 

.3223 

4 

3 

.4752“^ 

2 

2 

. 3299"^ 

4 

4 

.1570 

2 

3 

.5957'^ 

4 

7 

.995$"^ 

2 

4 

.1614 

4 

8 

.1105 

2 

5 

.1427 

4 

9 

,8203'^^ 

2 

6 

.8303"^ 

4 

10 

.4614"^ 

I 

I 
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Table  A-5 

MOLECULAR  CONSTANTS 

A 


State 


V 

B  (cm 
v'  ' 

D^(10“® 

F 

max 

Go  (v)  (cm 

0 

1.922 

6.0 

19220 

0.0 

1 

1.902 

6.0 

19020 

2174. 8 

2 

1.879 

6.0 

18790 

4317.0 

3 

1.861 

6.0 

18610 

6426.4 

4 

1.841 

6.0 

18410 

8502. 8 

5 

1.826 

6.0 

18260 

10545.8 

6 

1.808 

6.0 

18080 

12555. 0 

7 

1.781 

6.0 

17810 

14530. 7 

8 

1.766 

6.0 

17660 

16470. 6 

9 

1.740 

6.0 

17400 

18377. 6 

10 

1.724 

6.0 

17240 

20249. 7 

11 

1.703 

6.0 

17030 

22086.5 

12 

1.683 

6.0 

16830 

23887.6 

13 

1.663 

6.0 

16630 

25652.4 

14 

1.641 

6.0 

16410 

27380.5 

15 

1.620 

6.0 

16200 

29071.6 

16 

1.593 

6.0 

15930 

30725.1 

17 

1.572 

6.0 

15720 

32340.5 

18 

1.552 

6.0 

15520 

33017.3 

0 

2. 703 

6.0 

20730 

0.0 

1 

2.049 

6.0 

20490 

2371.5 

2 

2.025 

6.0 

20250 

4690. 3 

3 

2 . 002 

6.0 

20020 

6950. 7 

4 

1.968 

6.0 

19680 

9147.1 

S 

1.926 

6.0 

19260 

11269, 9 

6 

1.896 

6.0 

18960 

13310, 9 

7 

1.852 

6.0 

18520 

15262. 0 

8 

1.810 

6.0 

18106 

17100.2 

9 

1.762 

6.0 

17620 

18827. 1 

10 

1.710 

6.0 

17100 

20423.8 

11 

1.653 

6.0 

16530 

21903. 3 

12 

1.595 

6.0 

15950 

23275. 1 

13 

1.545 

6.0 

15450 

24551.4 

14 

1.494 

6.0 

14940 

25747. 7 

15 

1.452 

6.0 

14520 

26874.3 

16 

1.404 

6.0 

14040 

27941.4 

17 

1.355 

6.0 

13550 

23956.9 
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fable  A-6 


2-57-62-2 


I 

J 

I 

I 


I 

I 


FRANCK-CONDON  FACTORS  FOR  THE 
(FIRST  NEGATIVE)  SYSTEM  nJ 


rl 

V'-' 

q 

V 

V" 

q 

0 

0 

.  6509 

6 

4 

.  2673 

0 

1 

<  2588 

6 

5 

.  2830 

0 

2 

. 7016^1 

6 

6 

. 5331“1 

X 

0 

.  3014 

6 

8 

lOOR 

i 

1 

.  2226 

6 

9 

.1150 

1 

2 

.  2860 

6 

10 

. 7442*1 

1 

3 

.  1324 

6 

11 

. 3651*1 

1 

4 

*4273"‘l 

7 

4 

. 2789*1 

2 

6 

.4537^1 

7 

5 

.  3668 

2 

1 

.4060 

7 

6 

.2414 

2 

2 

. 5065‘1 

7 

7 

.7236*1 

2 

3 

.2290 

7 

9 

. 8084*1 

2 

4 

.1653 

7 

10 

.  1684 

2 

5 

. 7113*1 

7 

11 

. 8051*1 

3 

i 

.  1056 

7 

12 

.4411*1 

3 

2 

.  4137 

8 

5 

.  3419*1. 

3 

4 

.  1557 

8 

6 

.3401 

3 

5 

,  1706 

8 

7 

.  2083 

3 

6 

, 9451*1 

8 

8 

.  8468''^ 

3 

7 

. 3801"i 

8 

10 

, 5865*1 

4 

2 

.  1660 

8 

11 

. 9821*1 

4 

3 

.  3792 

8 

12 

. 8319*1 

4 

5 

, 9290*1 

8 

13 

. 5053*1 

4 

6 

.1569 

9 

6 

. 3890*1 

4 

7 

,1096 

9 

7 

.3686 

4 

8 

, 5236*1 

9 

8 

.1837 

5 

3 

,  2205 

9 

9 

, 9065*1 

5 

4 

,3310 

9 

11 

,4075*1 

5 

6 

,4815'-1 

9 

12 

. 8628*1 

5 

7 

.  1333 

9 

13 

. 8288=1 

5 

8 

.  1161 

9 

14 

.5544*1 

5 

9 

,6488*1 

10 

7 

.4153*1 

A-15 
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Table  A-6  (Continued) 


v' 

vt* 

q 

V 

V'* 

q 

10 

8 

.  3932 

12 

18 

. 4058*1 

10 

9 

.  1672 

13 

10 

. 3432*1 

10 

10 

.  9:141“1 

13 

11 

.  4496 

10 

13 

. 7393“^ 

13 

12 

.  1593 

10 

14 

.8018“i 

13 

13 

.7384*1 

10 

15 

. 5872-^ 

13 

16 

.4189*1 

10 

16 

. 3410“^ 

13 

17 

. 6404“^ 

11 

8 

. 4170“! 

13 

18 

. 5968*1 

11 

9 

.  4149 

14 

12 

.4614 

11 

10 

.  1579 

14 

13 

.  1697 

11 

11 

.  8‘818‘^ 

14 

14 

.6418*1 

11 

14 

.6208*1 

14 

17 

.3390*^ 

11 

15 

. 7574^1 

14 

18 

. 5774*1 

11 

16 

. 6041*1 

15 

13 

.4681 

12 

9 

. 3928"! 

15 

14 

.  1866 

12 

10 

.4338 

15 

15 

. 5355*1 

12 

11 

.  1554 

16 

14 

.4676 

12 

12 

. 8204*1 

16 

15 

.2104 

12 

15 

.  5131''^ 

16 

16 

.4236*1 

12 

16 

. 7018*1 

17 

15 

.4579 

12 

17 

.6066~1 

17 

16 

.2413 
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fable  A-7 


NO  MOLECULAR  CONSTANTS 


V 

B  (em"^) 

V  '  ' 

Y 

V 

Gl  (v)  (cm"^) 

Gg  (v)  (cm“ 

0 

1.6957 

5  X 

10"® 

73.24 

948.52 

948. 34 

1 

1.6779 

5  X 

16-6 

74.02 

2824.61 

2824.08 

2 

1.6601 

5  X 

10-6 

74.81 

4672.74 

4671.86 

3_ 

_ L  6423  __ 

_  S  X 

10-6 

75. 62 

_ ^2 

6491 , 70 

4 

1.6245 

5  X 

10"® 

76.46 

8285.13 

8283.56 

5 

1.6067 

5  X 

10“® 

77.  30 

10049. 37 

10047.44 

6 

1.5889 

5  X 

10“6 

78. 16 

11785.63 

11783.35 

7 

1.5711 

5  X 

10-6 

79.  05 

13493.91 

13491. 28 

8 

1.5533 

5  X 

10-6 

79.95 

15174.18 

15171.20 

9 

1. 5355 

5  X 

10-6 

80.  88 

16826.46 

16823.13 

16 

1.5177 

5  X 

10-6 

81.83 

18450.  73 

18447.06 

11 

1.4999 

5  X 

10-® 

82.80 

20047. 00 

20042.82 

12 

1.4821 

5  X 

1 

o 

83.  80 

21615.13 

21610.53 

13 

1.4643 

5  X 

10-6 

84.81 

23154. 94 

23149. 82 

14 

1.4465 

5  X 

10-6 

85.86 

24666.29 

24660. 64 

15 

1,4287 

S  X 

10-6 

86.93 

26148. 95 

26142. 79 

16 

1.4109 

5  X 

10-6 

88.02 

27602.67 

27595.91 
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Table  A-7  (Continued) 


State 

V 

®v 

D 

V 

G^iv) 

G2(V) 

B^n 

0 

1.118 

6  X  10“® 

28.6 

516.58 

517.30 

1 

1.105 

6  X  10^® 

29 

1638.66 

1540. 8 

y 

1.093 

_j6  X  10"®  _ 

31 

2546.4 

2550.0 

3 

1.081 

6  X  10-6 

32 

3540.4 

3545.4 

4 

1.068 

6  X  10"® 

29 

4521.2 

4527.7 

5 

1.056 

6  X  10"® 

36 

5489.4 

5497.3 

6 

1/041 

6  X  lO"® 

38 

6445.6 

6455. 1 

State  V 

D 

G(v) 

G^(v) 

®v 

V 

o 

A^S"-  0 

1.9870 

6  X  10 

1182.0 

0 

1 

1.9688 

6  X  10^® 

3523.4 

2341.4 

2 

1.9498 

6  10^® 

5833.4 

465i.4 

3 

1.9290 

6  X  10^® 

8110.2 

6928.2 

4 

1.9108 

6  X  10-® 

10352.® 

9170.7 

5 

1.8906 

6  X  10"® 

12562.4 

11380.4 

6 

1.8684 

6  X  10"® 

14740, 8 

13558.8 

7 

1.8486 

6  X  10'® 

16885,1 

15703.1 
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Table  A-8b 

FRANCK-CONDON  FACTORS  FOR  NO  GAMMA  SYSTEM  (A  ^2“^  -  X  ) 


v" 

q 

v' 

v" 

q 

0 

/1653 

4 

0 

.5073** 

1 

.2636 

1 

.2374 

2 

.2376 

3 

.  1263 

3 

.  1604 

5 

.5476“* 

4 

.9071^^ 

6 

. 8784“* 

10 

.5236“* 

0 

.  3291 

11 

.7116“* 

1 

.1049 

12 

.6971“* 

3 

.7214*^ 

13 

.5636** 

4 

.1349 

5 

.1339 

5 

1 

.1338 

6 

.  9858''^ 

2 

.  1976 

7 

. 6077*^ 

4 

.  1026 

5 

.5125** 

0 

.2906 

7 

. 6288** 

2 

.1545 

8 

.6791** 

3 

. 7523'^ 

13 

.6049“* 

6 

. 8848'* 

14 

.6148** 

7 

.1056 

8 

.8995~* 

6 

2 

.1989 

9 

.6292** 

3 

. 9697“* 

4 

. 5754'* 

0 

.1504 

6 

.8842“* 

1 

.1903 

9 

.6284’* 

4 

.1131 

10 

.5168'* 

5 

.5129** 

8 

. 6550** 

7 

2 

. 9467'* 

9 

.8549'* 

3 

.2113 

10 

.7933** 

5 

.  1058 

11 

.6037** 

7 

.7327“* 

11 

.5772'* 
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fable  A“9 


0.  MOLECULAR  CONST  ANT  S* 

A 


State 

V 

1  (cm 

D  (10  ^cm 

V  ' 

F 

max 

G  (v)  (em 

6 

1.438 

4.913 

9110,5 

0.0 

1 

1.422 

4. 825 

9011.8 

1556.4 

2 

1.406 

4. 737 

8913. 2 

3089. 1 

O 

X.otfU 

- CTtj±^,  O  — 

4698.3 

4 

1.375 

4.561 

8716.0 

6084.4 

5 

1. 359 

4.473 

8617.3 

7547.4 

6 

1.343 

4.385 

8518. 7 

8987.5 

7 

1.327 

4.297 

8420. 1 

10404. 9 

8 

1.311 

4. 209 

8321.4 

11799.7 

9 

1.296 

4. 121 

8222. 8 

13171. 8 

10 

1.280 

4, 033 

8124.2 

14521.4 

11 

1,264 

3;.  945 

8025.5 

15848. 3 

12 

1.248 

3. 857 

7926. 9 

17152.7 

13 

1.232 

3. 769 

7828.3 

18434,2 

14 

1.217 

3.681 

7729.6 

19693.0 

15 

1.201 

3.593 

7631.0 

20928.7 

16 

1, 185 

3.505 

7532.4 

22141.2 

17 

1, 169 

3,417 

7433. 8 

23330.3 

18 

1,154 

3.329 

7335.1 

24495,7 

19 

1,138 

3,241 

7236,5 

25637.2 

20 

1,122 

3,153 

7137-9 

26754.4 

21 

1.106 

3.065 

7039. 2 

27846,9 

22 

1,090 

2.977 

6940.6 

28914.4 

23 

1.074 

2.889 

6842.0 

29956.5 

24 

1, 059 

2. 801 

6743.3 

30972,6 

25 

1.043 

2.713 

6644,7 

31962.3 

26 

1,027 

2.625 

6546, 1 

32925,0 

27 

l.Oll 

2,537 

6447,4 

33860. 2 

28 

0. 9956 

2.449 

6348.8 

34767. 3 

29 

0,9798 

2.361 

6250.2 

35645. 7 

-1, 
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Table  A-10 

Og  SCHUMANN-RUNGE  FRANCK-GONDQN  FACTORS 


i! 

6 

0 

II! 

8 

9 

q(v'.v’') 

. 1253^^ 

»2802'‘^ 

£ 

1 

1 

tl 

13 

15 

1 

q(v',v") 

.2184’^ 

. 1833’^ 

U 

0 

10 

11 

.8667“^ 

1 

1 

17 

.  1137 

6 

12 

.1213 

1 

18 

.1301 

6 

13 

.  1465 

1 

19 

.1120 

0 

14 

.  1528 

2 

2 

.6472^'^ 

0 

15 

.  1379 

2 

3 

.4520''^ 

6 

16 

.1074 

2 

4 

.2200“^ 

0 

17 

.7210^^ 

2 

5 

.7886*^- 

0 

18 

*4159'^ 

2 

6 

.2138'^ 

0 

19 

.205r^ 

2 

7 

.4420’^ 

1 

2 

. 1294^® 

2 

8 

. 6899' - 

1 

3 

.  1022''^ 

2 

9 

.7819'“ 

1 

5 

.2407''^ 

2 

10 

.5804'- 

1 

6 

.7886’- 

2 

11 

.2000'^ 

1 

7 

.2049’- 

2 

13 

.2139'” 

1 

8 

.4262''- 

2 

14 

,6073’“ 

1 

9 

,7084’- 

2 

15 

.6959'” 

1 

10 

.9271'” 

2 

16 

,3556"” 

1 

11 

.9195’^ 

2 

18 

. 1536’^ 

1 

12 

.6286’^ 

2 

19 

.7004'^ 

A<^23 
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2-57-62-2 


V' 

q  (v',v") 

V 

xl! 

q(v'.v") 

3 

0 

. 1043“^ 

4 

6 

*5306"^ 

3 

1 

. 2243-^ 

4 

7 

.6118"^ 

3 

2 

.2248"^ 

4 

8 

.3937"^ 

3 

3 

. 1386“^ 

4 

11 

.3550"^ 

3 

4 

. 5844"^ 

4 

12 

. 4551“^ 

3 

5 

.1767“^ 

4 

13 

. 1519"^ 

3 

6 

,3896“^ 

4 

15 

.3357"^ 

3 

7 

.6189“^ 

4 

16 

.4976“^ 

3 

8 

.6716“^ 

4 

17 

.1687“^ 

3 

9 

.4256*^ 

4 

19 

.4176*^ 

3 

12 

.3544"^ 

5 

1 

.  17  ir^ 

3 

13 

.5616'’^ 

5 

2 

. 1374"^ 

3 

14 

.3187^" 

5 

3 

*6558"^ 

3 

16 

. 1776"^ 

5 

4 

.2038"^ 

3 

17 

.5938^^ 

5 

5 

,4224"^ 

3 

18 

.5971''- 

5 

6 

. 5654"^ 

3 

19 

.  mo’^ 

5 

7 

.4268*“ 

4 

0 

.3479'"® 

5 

8 

. 1022"" 

4 

1 

,6790'^ 

5 

10 

.2965"^ 

4 

2 

, 6084" " 

5 

11 

.4003"^ 

4 

3 

.3310"^ 

5 

12 

, 1149"^ 

4 

4 

. 1202"^ 

5 

14 

.3425"“ 

4 

5 

.3033"^ 

5 

15 

.36U"^ 
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T^le  A-10  (Continued) 


t 

XH 

q  (v',v") 

x! 

xl 

q  (v’.v") 

5 

16 

«3720"^ 

7 

4 

.3759"^ 

5 

17 

. 1360"^ 

7 

5 

.4900"^ 

5 

18 

.4761^^ 

7 

6 

. 3094"^ 

_-I 

5 

19 

.2836“" 

7 

9 

.3321 

6 

1 

.3736“^ 

7 

10 

.2094“^ 

6 

2 

.2683"^ 

7 

12 

.2098“^ 

6 

3 

.1121“^ 

7 

13 

.3004"^ 

6 

4 

.2956“^^ 

7 

15 

. 1305'^ 

6 

5 

.4930“^ 

7 

16 

.3437-^ 

6 

6 

.4767^^ 

7 

17 

.8005“^ 

6 

7 

a874"^ 

7 

18 

.8667"’^ 

6 

9 

.aoas'^^ 

7 

19 

.3852"^ 

6 

10 

,3754^^ 

8 

2 

.7343“^ 

6 

11 

. 1415^^ 

8 

3 

.2323"- 

6 

13 

.2919'- 

8 

4 

.4259"- 

6 

14 

.3050'^ 

8 

5 

.4158*- 

6 

16 

.1677”^ 

8 

6 

. 1390"- 

6 

17 

.4902'^ 

8 

8 

.2307"- 

6 

18 

, 1169"- 

8 

9 

.2829’- 

6 

19 

.6046“^ 

8 

11 

. 1033"^ 

7 

1 

.7262"^ 

8 

12 

.2975"“ 

7 

2 

.4660"- 

8 

15 

, 3032’- 

7 

3 

. 1698"^ 

8 

16 

.1014’- 

A-25 
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Table  A- 10  (Continued) 


v' 

y" 

q(v*.v") 

I' 

xl! 

q(v’,v") 

8 

17 

.5726"^ 

10 

11 

-  -  -“1 

.2386 

8 

18 

.3281“^ 

10 

13 

.1271"^ 

8 

19 

. 1098*^ 

10 

14 

.2247’^ 

10 

ic 

_ ,  1795’^ 

9 

— - 2 — 

•  XUOu 

9 

3 

.2910*^ 

10 

17 

. 1902’^ 

9 

4 

.4341“^ 

11 

3 

.3645’^ 

9 

5 

.2967"^ 

11 

4 

.3327’^ 

9 

8 

.2931"^ 

11 

5 

.6897“^ 

9 

9 

, 1313“^ 

11 

7 

.2481“^ 

9 

11 

.2422"^ 

11 

8 

. 1333’^ 

9 

12 

. 1661“^ 

11 

10 

.2122’^ 

9 

14 

.2382^^ 

11 

11 

.1208"^ 

9 

15 

. 1607^^ 

11 

13 

.2231*^ 

9 

16 

. 1294'^ 

11 

16 

.2304"- 

9 

17 

*2692^^ 

12 

3 

.3695’^ 

9 

18 

, 1339”^ 

12 

4 

.2477"^ 

10 

2 

. 1439"^ 

12 

6 

. 1126’- 

10 

3 

.3372’^ 

12 

7 

.2375’- 

10 

4 

.4000’“ 

12 

10 

.2114’^ 

10 

5 

.  1701"“ 

12 

12 

. 1175’- 

10 

7 

, 1903’“ 

12 

13 

.1758”- 

10 

8 

.2433"“ 

12 

IS 

.1868’- 

10 

10 

.1193"^ 

13 

3 

.3526’^ 
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Table  A*  10  (Continued) 


v' 

q(v',v") 

v' 

q(v'.v") 

13 

4 

*  1616' ^ 

16 

5 

.  IMO"^ 

13 

6 

.1764“^ 

16 

6 

. 1694“^ 

13 

7 

. 1750“^ 

16 

8 

4 1274"^ 

13 

9 

. 1582“^ 

16 

9 

. 1062“^ 

13 

10 

.1339“^ 

16 

11 

.1597“^ 

13 

12 

. 1878  ^ 

16 

14 

. 1066"^ 

13 

15 

*  1854“^ 

17 

5 

. 1484“^ 

14 

4 

i  3176"’^ 

17 

6 

. 1199“^ 

14 

4 

.8805"^ 

17 

8 

. 1536“ - 

14 

6 

17 

11 

. 1330“^ 

14 

7 

,9600"^ 

17 

13 

. 1393“^ 

14 

9 

.  1898''^ 

18 

5 

.166r- 

14 

12 

. 1685"^ 

IS 

8 

.149r^ 

14 

14 

,1528’- 

18 

13 

. 1328““ 

15 

6 

.2040^^ 

19 

5 

. 1661“^ 

15 

9 

. 1649"^ 

19 

8 

. 1210““ 

15 

11 

. 1358'^ 

20 

5 

. 1515"^ 

15 

14 

.1673'^ 
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AppeMix  1 

SUPPLEMENTARY  NUCLEAR  WEAPON  EFFECTS  STUDY 

B.  1  INTRODUCTION 

The  proper  theoretieal  description  of  the  fireball  resulting  from  a  nuGlear  detonation 
requires  that  the  radiative  transport  of  energy  be  ineluded  in  the  appropriate  hydros- 
dynamieal  description  of  the  fireball  development.  When  this  is  attempted^  it  is 
quickly  found  that  numerical  calculations  are  required  to  solve  the  equations  involved. 
The  machine  codes  for  such  numerical  computations  are  commonly  referred  to  as 
coupled  radiation ‘-hydrodynamics  codes. 

Brode  (Refi  23)  j  of  the  Rand  Corporation,  has  used  such  codes  for  a  variety  of  fireball 
computationSi  His  results  are  valid  at  earlier  times  when  the  fireball  is  essentially 
Optically  thick  and  when  the  energy  lost  is  small  eompared  to  tliC  internal  energy  of 
the  system.  In  terms  of  low  altitude  fireball  phenomenology,  his  results  are  valid  at 
times  earlier  than  the  time  to  light  minimum.  It  is  to  be  noted  that  the  fireball  is  still 
optically  thick  at  this  time,  and  only  a  small  fraction  of  the  weapon  yield  has  escaped 
the  system  via  thermal  radiation.  It  should  further  be  noted  that  the  time  of  minimum 
corresponds  approximately  to  the  time  at  wMch  the  shock  wave  departs  from  the  radia-^ 
ting  region  of  the  fireball,  and  the  Brode  results  should  thus  give  an  adequate  descrip.- 
tion  of  the  advancing  shock  front. 

Proper  theoretical  description  of  the  final  pulse  of  thermal  radiation,  which  involves 
as  much  as  half  of  the  weapon  yield,  requires  that  the  radiation-hydrodynamics  code 
be  capable  of  treating  radiative  transport  in  botii  optically  tijick  and  optically  Uiin  media 
in  conjunction  with  the  appropriate  hydrodynamics.  Such  a  code  has  been  develt^ed 
and  is  described  below. 
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B.  2  GENERAL  DESGRIPTION 


The  code  developed  for  late  tiffle  firehail  computations  is  a  generalization  of  the 
Riehtmyer=von  Neumann  artificial  viscosity  technique  for  a  spherical  Lagrangian 
system  in  which  the  divergence  of  the  radiative  flux  appears  in  the  hj^erbolic  partial 
differential  equation  which  expresses  the  conservation  of  energy  for  a  given  mass  zone^ 
The  basic  dilferential  equations  are  those  given  by  Ricbtmyer  (Ref^  24) 
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modified  to  include  this  additional  radiative  term.  Without  this  term,  the  solution  of 
the  equations  is  a  standard  initial  value  problem  which  can  be  integrated  explicitly. 


Such  solutions  are  well  known;  only  those  modifications  required  to  include  the  radia¬ 


tive  transport  will  be  discussed  here. 


Formulation  of  the  Radiative  Energy  Transport 

To  understmid  the  basic  problem  of  radiative  transport, consider  the  standard  equation 
of  radiative  transfer  for  a  medium  in  local  thermodynmnic  equilibrium  and  for  a  plane 
parallel  geometry: 

CPS  e  ^  |b^  jT(x)|  -  I^(x  cos  d)|  (B,4) 

where  the  quantities  are  defined  in  the  standard  termhiology  as  given  by  Chandrasekhar 
(Ref.  25). 
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This  equation  represents  the  change  in  the  specific  intensity  of  the  radiation  along  a 
ray  having  a  given  direction  in  spacei  Since  an  infinite  nuniber  of  directions  for  rays 
can  be  assigned  to  any  unit  volume  imder  consideration,,  it  is  clear  that  Eq«  (B^4) 
represents  an  infinite  set  of  differential  equations  which  are  functions  of  position, 
direction  in  space,  and  radiation  frequency. 

The  solution  for  this  set  of  equations  has  been  obtained  in  closed  forna  by  Chandrasekhar 
-^efi  25  ,  55  f^  based  on  a  method  4n  whieh^the  radiation  field  is  represented  by  — - 

2N  streams  where,  in  the  limit,  N  approaches  infinity.  The  convergence  of  Ihis 
solution  is  shown  to  be  critically  dependent  on  the  proper  choice  of  radiation  streams 
used  in  the  first  approximation,  i.  e. ,  N  =  1,  which  corresponds  to  the  twO-stream 
method  of  A.  Schuster  and  K.  Schwarzschfld  (1905).  Suitable  choice  of  the  directions 
of  the  two  rays  used  to  represent  the  radiation  field  results  in  a  formulation  of  the 
gross  energy  transport,  whose  accuracy  is  considerably  greater  than  the  acGUraey 
with  which  we  know  the  appropriate  absorption  coefficients .  The  accuracy  Of  this 
{qtproach  improves  when  it  is  applied  to  thin  isothermal  slabs  as  contrasted  with  the 
usual  application  to  an  atmosphere  of  essentially  infinite  extent. 

The  basic  physics  of  the  Schuster ^chwarzschild  method  has  been  carried  over  to 
spherical  geometry.  In  this  application,  the  inclination  of  the  chosen  radiation  streams 
changes  with  the  relative  radius  of  curvature  of  the  spherical  zones.  When  the  radius 
of  curvature  becomes  large  compared  to  the  zone  thickness,  the  appropriate  inclination 
is  that  which  is  applicable  for  plane  ^parallel  geometry.  Indeed,  in  this  limit,  the 
exact  solution  for  the  flux  in  terms  of  exponential  integrals  may  be  employed  and  no 
approximation  is  required.  However,  in  the  present  application,  other  conditions  of 
the  problem  made  it  desirable  to  maintain  the  two^stream  approximation,  the  choice 
of  inclination  being  determined  by  the  condition  that  the  difference  in  these  two  solutions 
be  held  to  a  minimum.  The  proper  choice  of  inclination  near  the  center  of  the  sphere 
has  been  temporarily  postponed  by  choosing  the  central  zones  sueh  that  they  are  optic* 
ally  thick  and  by  requiring  that  the  two*stream  solution  transform  into  the  diffusitm 
approximation  vdiich  is  sqiplicable  under  these  conditions. 


B*3 


LOCKHEED  MISSILES  tk  SPACE  COMPANY 


2.57-«62--2 


The  actual  formulation  of  the  flux  is  then  carried  out  using  the  formal  integral  solution 
(Ref*  25)  of  the  equation  of  transfer  for  the  specific  intensity  along  the  representative 
rays*  The  source  function  is  assumed  to  be  that  which  is  appropriate  under  conditions 
of  local  thermodynamic  equilibrium  and  is  expanded  in  ter  ms  of  ibe  temperatures  — 
and  derivatives  of  the  temperatures  —  at  the  zone  boundaries.  While  the  code  has  been 
written  to  allow  the  spectrum  to  be  divided  into  any  desired  number  of  frequency  bands, 
the  initial  Computations  have  been  Garried  out  using  the  Rosseland  mean  absorption 
coefficients  for  the  purposes  to  be  enumerated  later. 

The  numerical  solution  is  carried  out  by  explicit  integration  of  the  basic  equations, 
except  for  the  energy  equation,  using  die  well^lmown  Richtmyer-von  Neumann  method 
(Ref.  26).  The  integration  of  the  energy  equation  cannot  be  carried  out  ejqilieitly  since 
the  radiative  flux  term  for  any  particular  zone  is  dependent  upon  the  temperatures  of 
the  neighboring  zones.  Since  the  density  is  known  at  this  stage  of  the  solution,  it  is 
convenient  to  ej^ress  the  energy  equation  in  terms  of  zone  temperatures  alone.  The 
solution  is  then  carried  out  by  assuming  an  ^proximate  temperature  for  each  zone 
and  using  a  Newton^Raphson  iteration  scheme  to  adjust  the  zone  temperatures  until 
changes  in  temperatures  on  subsequent  iterations  are  less  than  a  pre -determined 
amount. 

This  method  of  iteration  involves  the  solution  of  N  linear  equations  in  N  unknowns, 
where  N  is  the  number  of  zones  and  is  usually  of  the  order  of  a  hundred.  Such  a 
system  is  in  principle  easily  soluble  by  computer  methods,  but  the  procedure  is  time 
consuming.  In  practice  only  a  few  zones  in  each  direction  measurably  influence  the 
temperature  iteration  for  a  particular  zone,  which  results  in  a  simplification.  In  a 
first  solution  of  the  problem  a  five -zone  centered  system  has  been  used  which  results 
in  a  mateix  having  five  terms  across  the  diagonal'  After  suitable  scaling,  such  a 
system  is  rtmidly  soluble  by  direct  elimination  and  back  substitution. 

As  part  of  the  check  out  procedure,  the  equation  of  state  and  formulation  of  the  mean 
absolution  eoeffieient  given  by  Erode  (Ref.  23)  were  used  so  as  to  reproduce  bis  results 
over  an  apprqpriate  period  of  fireball  history.  Using  100  zones,  and  3  iterations  per 
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hydrddynamieal  cycle,  the  code  proceeds  successfully  about  200  time  steps  per  hour. 
With  minor  modifications,  the  code  has  also  been  applied  to  a  problem  in  thermal 
ablation  in  fireballs  and  to  problems  in  stellar  stability. 

Interpolation  formulae  are  currently  being  written  for  mean  absorption  coefficients 
for  chosen  spectral  bands  appropriate  to  late  fireball  development. 

Initial  Energy  Deposition  in  Variable  Density  Atmosphere 

When  a  nuclear  device  is  detonated  in  the  atmosphere,  a  large  percentage  of  the  eoer^ 
released  is  radiated  away  and  is  deposited  in  the  surrounding  atmosphere.  Since  the 
air  density  exhibits  a  large  variation  with  altitude ^  it  follows  that  the  scale  of  distances 
over  which  this  energy  is  absorbed  by  the  air  will  be  dpendent  upon  the  burst  altitude. 
Consequently  the  energy  density  also  exhibits  a  dependence  on  altitude. 


The  sequence  of  radiative  and  hydrodjmamic  events  which  follows  this  initial  phase  is 
dependent  upon  the  radial  energy  distribution.  Thus,  one  would  eJqpect  the  scaling 
laws,  for  various  weapons  effects  as  a  function  of  altitude,  to  also  exhibit  a  dependence 
on  the  initial  energy  deposit.  It  is  therefore  of  interest  to  perform  computations  of 
this  distribution  as  a  function  of  weapon  yield  and  altitude. 

Because  of  the  large  number  of  cases  of  interest,  a  machine  code  was  written  to  per¬ 
form  the  i^ropriate  computations.  The  basic  equations  involved  are 
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=  KiPtP)  (®*7) 

T  =  T(E,p)  (B.8) 


p  -  air  density  as  a  funetion  of  altitude  A 
-  absorption  coefficient  of  air 
W  -  total  energy  of  freqaency  radiated  by  the  source. 

Equation  (Bi  6)  was  evaluated  by  using  an  empirical  fit  to  the  ARDC  Standard  Atmosphere 
Olef.  27). 


The  absorption  coefficients  used  were  tiiose  of  Gilrnore  (Ref.  28)  for  cold  air.  Stripping 
was  neglected  since  it  appears  to  have  only  a  small  effect  upon  the  later  development. 


In  coding  the  problem,  the  radiative  output  of  the  source  was  divided  into  ll  separate 
spectral  regions,  and  a  corresponding  average  opacity  was  chosen  for  each  of  these 
regions.  A  minor  modification  of  the  code  allows  a  larger  number  of  spectral  regions 
to  be  used  when  desired.  The  code  carries  out  radial  computations  for  the  downward 
direction  and  at  directions  in  space  differing  by  10  deg^  or  any  multiple  of  10  deg,  as 
may  be  desired.  The  printed  output  for  each  direction  in  space  lists  the  radius  from 
the  point  of  detonation  for  35  different  points  and  the  air  density,  altitude,  and  energy 
density  for  each  point. 

Temperatures  were  determined  using  the  tabulated  equation -of -state  data  given  by 
Gilmore  (Ref.  16),  and  Hilsenrath  and  Beckett  (Ref.  29).  Initisd  runs  have  shown  that 
a  better  representation  of  the  equation  of  state  for  air  is  desirable. 
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The  initial  funs  inGluded  a  vafiety  of  yields  and  altitudes,  including  eases  applicable 
to  pfevious  field  test  situations  and  cases  of  intefest  for  future  planning.  Considerable 
success  was  obtained  in  predicting  a  nutnber  of  the  effects  observed  in  earlier  field 
tests  i  but  a  better  handling  of  the  equation  of  state  is  indicated  before  the  results  are 
presented  in  scaling  law  form. 
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